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ABSTRACT 


This thesis addresses the development of a physically based modeling simulator 
for a ship in a 3-D virtual environment to be used in naval tactical training systems. The 
objective is to develop a computer simulation program in which physical models are im¬ 
plemented in order to achieve a realistic representation of a ship in a virtual environment 
considering its physical features in the presence of environment conditions including 
waves, ocean current, wind, fog and day/night issues. The simulator was developed by in¬ 
tegrating five marine models with a virtual ocean environment created with a visual 
simulation builder tool. The marine models include a maneuvering model, a wave model, 
a wind model and an ocean current model. The numerical results from another complex 
wave model were also combined using linear interpolation to increase the realism level of 
the simulator. The result of this thesis shows that the integration of multiple models from 
different sources is a feasible approach to meet the application requirements. The result 
also indicates that the use of the interpolation technique to take advantage of complex 
models yields a simulator with acceptable level of realism while imposing very low com¬ 
putational load in the application program. 
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EXECUTIVE SUMMARY 


This thesis addresses the development of a physically based modeling simulator 
for a generic ship in an open water 3-D virtual environment to be used primarily in naval 
tactical training systems. The study includes marine dynamic models and their integration 
for the purpose of developing a simulator using virtual reality techniques. The main ob¬ 
jective is to develop a computer simulation program in which physical models are im¬ 
plemented in order to achieve a realistic representation of a ship in a virtual environment 
considering physical features of the ship such as engine, propeller and skull, in the pres¬ 
ence of environment conditions including waves, ocean current, wind, fog and day/night 
issues. 

The approach taken in this thesis was to integrate five marine models with a vir¬ 
tual ocean environment created with the support of a visual simulation builder tool called 
Vega from MultiGen-Paradigm Inc. Vega is an easy-to-use development tool that incor¬ 
porates modules for marine and special effects. The object of the simulation is a high¬ 
speed container ship. The implemented models include a 4-DOF (degrees of freedom) 
maneuvering dynamic model in surge, sway, roll and yaw, a 1-DOF wave model that 
adds random behaviors to the roll motion, and a 3-DOF wind model in surge, sway and 
yaw. A complex 6-DOF wave model was also combined with the other models to add the 
pitch and heave motions. This complex wave model runs as a separate application and 
generates numerical tables that were interpolated in real-time to obtain the six motion 
outputs during the simulation. The fifth model implemented in this study is a 2-DOF 
ocean current model in surge and sway motions. Some features were added to the simula¬ 
tion program to facilitate the observation of the environment forces and to verify the ship 
performance. They are a virtual grid on the ocean, a compass, arrows that represent the 
environment force vectors and a stopwatch. 

The principal objective of this thesis was reached. The simulation program 
worked well and this study showed that the integration of multiple models from different 
sources is feasible and produces good results, meeting the application requirements. The 
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use of the interpolation teehnique to take advantage of eomplex models worked well, 
produeing very low eomputational load in the applieation program with very aeeeptable 
results. The use of four dynamie models with relative low eomplexity produeed suffi- 
eiently good results, saving proeessing time for other tasks when running on eomplex 
simulation systems. The Vega system showed to be very useful to build eomplex virtual 
environment applieations in a short time. 



I. INTRODUCTION 


A. BACKGROUND 

As a result of the rapid and cheaper development of computer graphics and hard¬ 
ware, the use of virtual environments (VE) simulators has become a necessity. In military 
applications, particularly, it represents a feasible approach to experiment with new tactics 
and weapons, in order to interconnect training units worldwide affordably. 

This thesis addresses the preliminary development of a physically based modeling 
simulator for a generic ship in an open water 3-D virtual environment to be used primar¬ 
ily in naval tactical training systems. It could also be useful in many other maritime simu¬ 
lation areas, such as naval damage control simulators, marine vessel simulators for stabil¬ 
ity analysis and training systems for entering and leaving port. The complexity of the 
models to be used will depend on the intended applications. 

From the point of view of control systems, this thesis presents ship models with 
four and six degrees of freedom, environment models and interaction between the mod¬ 
els. 

When considering naval tactical simulators, the main goal is to take advantage of 
several features of VE techniques to increase the functionality of a naval tactical simula¬ 
tor. It is not only a matter of obtaining better visual presentation using many texture tech¬ 
niques and viewpoints possibilities. The point is that a simulator using VE must have 
more functionality than a simulator not using VE. In other words, in a VE tactical simula¬ 
tor, the user must have more resources to obtain useful information from the operation 
scenery to help him make decisions than when utilizing the old fashion simulators. 

For example, if a trainee officer, using the console/cage set as a helicopter, could 
see the image of a maneuvering ship in the ocean from the virtual cockpit, then the 
trainee officer would be able to infer which type of evasive movement the ship is under¬ 
taking by simply observing the ship’s wake. This kind of information is what is available 
to the user in a VE simulator that would allow the user to feel immersed in the operation. 
In this same sense, a virtual representation of the sea state, weather conditions, the ship 
rocking, and so forth, would add a new condition to the exercise. The attitudes taken by 
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an officer during good weather and flat sea scenery certainly would be different from the 
ones taken during a foggy situation or in a hard sea state condition where the visual clues 
are bad. Further, the pitch and roll movements of the ship in a stormy sea state condition 
may change some reeeption features of sensors, such as radars, thereby ereating effeets 
that are very diffieult to handle, such as the Sea Clutter, the Image Effeet, and the Sea 
Skimmer, for example. It may also impose eonstraints in some armament usage. 

As ean be seen, the functionality and the feeling of immersion would be inereased 
due to the intrinsic characteristics of a VE tactical simulator. 

Naturally, all sensor and weapon models used in the taetieal simulator would have 
to be concerned with this new input information, and their eomplexity would eertainly in- 
erease to support the roeking ship scenario. Therefore, the eomputational load of the sys¬ 
tem would also increase. The tradeoff between the model’s complexity and processing 
time will be discussed later in the following ehapter. 

B, PROBLEM STATEMENT 

The researeh question of this thesis is how to simulate physieal behaviors of a 
ship in an open-water virtual environment with respect to its structural characteristics and 
the environment conditions. 

This thesis addresses marine dynamie models and their integration for the purpose 
of developing a VE tactical simulator. More specifieally, the objeetive is to develop a 
computer simulation program where physical models are implemented in order to achieve 
a realistic representation of a ship in a virtual environment considering its physical fea¬ 
tures such as the engine, propeller and skull, in the presence of environment conditions 
including waves, ocean currents and wind, fog and day/night issues. 

The models used in this study are obtained from the literature. Some adaptations 
are neeessary to work with different models from different sources for the same ship. 
This is a result of the laek of complete information about the ship that was chosen for this 
simulation from any single source. 


2 



C. THESIS OUTLINE 

Following this introduction, Chapter II describes how the objectives of this simu¬ 
lation are obtained in terms of what is more relevant to model in what kind of application. 
Some considerations about the tradeoff between the complexity of the models and real¬ 
ism complete this chapter. Chapter III offers an overview of the simulation system pre¬ 
senting its main features and the system modules. Chapter IV describes the Virtual Envi¬ 
ronment Module. The marine module and the marine special effects module provided by 
the visual simulation system are also presented. Chapter V explains the Ship’s Motion 
Module. A theoretical background in hydrodynamics is first discussed to introduce some 
notation and concepts that will make it easier to understand the models. Then, each model 
is described, starting with the Maneuvering Model. Finally, the environment disturbance 
models including the Wave Model I and II, the Wind Model and the Ocean Current 
Model are presented. Chapter VI describes the implementation of the simulation program 
and explains all the principal modules. Chapter VII presents the method to operate the 
simulation program and Chapter VIII finalizes presenting the conclusions and recom¬ 
mendations for future work. 
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II. OBJECTIVES OF SIMULATION SYSTEM 


A, WHAT TO MODEL? 

There are many parameters to be eonsidered in developing a ship simulation sys¬ 
tem. It is not possible to consider all parameters in this thesis work. A natural question is: 
what would the users and experts in tactical training consider important modeling pa¬ 
rameters in terms of physical characteristics of the ships and the environment that could 
make the training more realistic and immersive and, consequently, would increase the 
functionality of the simulator? 

In order to answer the initial question and determine the most important modeling 
parameters in terms of physical characteristics of the ships and the environment, the fol¬ 
lowing approach was taken. First, some interviews were conducted to ascertain the main 
issues that could be explored to take advantage of the VE usage in the ship simulation. 
These issues are related to the modeling of behaviors of a ship under various weather/sea 
state conditions. 

Interview techniques are useful for identifying possible areas for more detailed 
analysis, and interviews are easy to conduct and direct. The unstructured interview can 
generate interesting points and statistical analysis can be run on the user’s answers [1]. 
The data collected provide information about general rules and principles and is faster 
than observational techniques. Interview techniques are useful for investigating events 
that occur infrequently. In addition, the interviews can be recorded for future analysis. 

1. Survey 

Based on the results of the interviews, a survey was created to obtain the opinions 
from a considerable number of users in order to rank modeling parameters in terms of 
relevance. The subjects of this survey were surface warfare officers from the Brazilian 
Navy and U.S. Navy, and some people involved with naval tactical training in Brazil. 
Their opinions were collected and analyzed in order to prioritize the research efforts of 
this thesis. 
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The survey used in this study appears in Appendix A, exactly as presented to the 
subjects. 

2. Results and Analysis 

The opinions of fifty people (50% from Brazil and 50% from the United States) 
were taken and analyzed. The analysis method used was the simple summation of the 
rank values of each option. Since the most relevant option receives the rank value 1 
(one), the option with the lowest summation value will be the most relevant. The results 
are shown in Table 1. 


Relevance 

Option 

Total 

1 

acceleration/deceleration rates 

92 

2 

turning rate 

93 

3 

visibility issues (day/night, fog) 

193 

4 

influence of sea state (ship “rocking”) 

195 

5 

influence of the wind 

201 

6 

influence of currents 

207 

7 

weight distribution (tilt/trim) 

284 

8 

wake of the ship 

296 


Table 1. Survey Results. 


The results demonstrate that users and experts established three distinct categories 
of relevance. The two first options had rank values considerably lower than the others. 
Thus, a realistic ship movement is basic in this kind of simulator. The second group in¬ 
cludes the influence of the environment on the ship behavior. The third group, with the 
least relevant values, is the weight distribution and the wake of the ship. 

The option “Other” was used with several secondary suggestions. Most were re¬ 
lated to the modeling of engines, such as the treatment of ship damages that could con¬ 
straint the performance of the ship in a warfare operation. This feature will not be con¬ 
templated in this study but it is an interesting and important area for future implementa¬ 
tions. 
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B, MODEL COMPLEXITY TRADEOFF 

Another important issue is the complexity of the models. It is expected that com¬ 
plex models are more expensive in terms of processing time. An important aspect to bear 
in mind is that models do not have to include all the details of the real world. Models are 
evaluated by their ability to provide reasonable answers. They do not have to be perfect 
representations of the real world but they need to reasonably represent those aspects of 
the real situation that greatly influence the measurements of interest. In this kind of appli¬ 
cation, there are visual and analytic outputs to evaluate the effectiveness of the simulation 
system. Therefore, the models need to be simple enough to create a good visual represen¬ 
tation contributing to an increase in the sense of immersion and be able to produce out¬ 
puts that interact well with other simulation models. 

In the process of system implementation, it is necessary to adjust the complexity 
of some algorithms according to the margin of spare time of processing of the system. 
This kind of problem is very common when the system relates to networked virtual envi¬ 
ronment applications [2]. Many times the system must be “tuned” to obtain a better time 
distribution among rendering, computational calculations and communication processes. 

In order to achieve good results combining simplicity and realism, an interpola¬ 
tion technique was tried in this simulation, specifically, in the ship’s wave model. This 
technique allows the use of the results of a very complex wave model in a very cheap 
computational way, by interpolating the desired parameters from tables previously gener¬ 
ated. 

This chapter described how the main objectives of this thesis were achieved based 
on the opinions of expert people. The next chapter provides an overview of the simula¬ 
tion system by describing how it expects to achieve the desired objectives. The principal 
features of the simulator are then presented and its component modules are discriminated. 
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III. SIMULATION SYSTEM OVERVIEW 


A, MAIN FEATURES 

Based on the results of the survey, the objective of this thesis is to build a pre¬ 
liminary simulation block to possibly observe the behavior of a ship under some envi¬ 
ronmental conditions. The results of this study would be used as a part of a larger simula¬ 
tion system where other capabilities would be included. Therefore, the main feature of 
this simulation system is to provide the user ship maneuvering capabilities and to allow 
changes in the environment conditions. A simulation is presented in a virtual scene in a 
graphic terminal. 

The maneuvering set includes the acceleration/deceleration of ship (engine order) 
and turnings (rudder order). The environmental settings make it possible to define the 
amount of daylight in the scenario, fog intensity, the sea state level and the wave-heading 
angle, the direction and intensity of the wind, and the direction and intensity of the ocean 
current. Weight distribution issues are not included in the scope of this thesis. 

The simulation system provides real-time information about the ship’s position 
and orientation, current engine state, shaft speed and rudder angle, direction and intensity 
of the wind and ocean current vectors, and the wave-heading angle. This information may 
be textual and/or graphical. 

It is possible to vary the observer’s position. The options include the following 
points of view: perspective view, lateral view, fore view, astern view, above view and on- 
the-bridge view. 

B, SYSTEM MODULES 

The simulation system is divided into modules. In general, it is very unusual in 
this type of simulation to have a single block including all the modeling equations. Espe¬ 
cially in this study, where the principal objective was not to model but to use and to inte¬ 
grate existing maneuvering and environmental disturbance models, the modularity was a 
necessity. Based on this fact, the simulator is composed of two main modules. 
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The Virtual Environment module holds three sub-modules. The User’s In¬ 
put/Output is in charge of accepting commands that set internal parameters associated 
with the ship maneuvering and environment related models, and displaying textual in¬ 
formation on the screen. The Graphics sub-module is responsible for rendering the virtual 
reality scene. The Ocean Model and Special Effects sub-module contains the set of equa¬ 
tions that represent the ocean dynamics, ship’s wake and bow-wave, field vectors such as 
wind and ocean current, and some special effects such as ship’s smoke. This module also 
treats daylight and fog issues. 

The Ship Motion module holds four sub-modules. Each sub-module plays a spe¬ 
cific role in the dynamics that contribute to representing the expected ship behavior. The 
Maneuvering Model sub-module treats the acceleration/deceleration and the turnings in¬ 
puts. The Wave Model I and II, the Wind Model and the Ocean Current Model sub- 
modules provide the motions produced by the sea state conditions, wind and ocean cur¬ 
rents respectively. Eigure 1 shows the block diagram of the system. 


Simulation 

System 


>■ 


t 



Virtual Environment Module 


Ocean Model & Special Effects 


User's 

I/O 


Graphics 


Maneuvering Model 


Wave Models I & II 


Wind Model 


Ocean Current Model 


Ship Motion Module 


Eigure 1. System Modules. 
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C. VALIDATION OF THE MODELS 

Each model of the Ship Motion module is subjeeted to a validation proeess. The 
validation proeess has two basic approaches, the analytie and the visual tests. The ana- 
lytie test is based on measurements of variables of interest that are relevant to a specific 
test. The measured variables are plotted to make their analysis easier. All measurements 
were made with the aetual simulation program running. A trigger eommand was imple¬ 
mented in the user’s interfaee to start the reeording of pre-defined variables of interest in 
a textual log file at any time. This log file is in a matrix format in whieh eaeh row ele¬ 
ment is related to a measured variable on the simulation test. The first row element is 
usually reserved for the time variable. Therefore, at eaeh iteration step, a new row is ap¬ 
pended, eontaining the eurrent time value and the respective measurements. When the 
trigger eommand is hit a seeond time, the log file is elosed. Then, the log file is imported 
to the Matlab workspaee and the plots are ereated. The analyses are made through those 
plots. 

The visual test is made by direet observation of the ship behavior in the virtual 
environment. Basically, this test is based on the eommon sense about how a ship is sup¬ 
posed to behave under eertain eonditions. 

After the diserimination of the modules, the following step eonsists of detailing 
the Virtual Environment module. The next ehapter deseribes the main eapabilities of the 
visual simulation system used in this study. 
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IV. THE VIRTUAL ENVIRONMENT MODULE 


There are many tools to help in the development of virtual reality systems due to 
the high eomplexity of the applieations. To start the rendering using only the original 
OpenGL functions is practically impossible. The tool used in this thesis is called Vega, 
developed by MultiGen-Paradigm Inc. Vega has been used in many projects in the Mod¬ 
eling, Virtual Environments and Simulation Department (MOVES) at the Naval Post¬ 
graduate School, supporting different virtual reality applications. 

A, THE VEGA SYSTEM 

1. General Description 

Vega is a high performance visual simulation system [4]. It consists primarily of 
an executable application, several C callable libraries and a graphical user interface. It is 
modular and extensible in design, and provides users a simple and direct way of interfac¬ 
ing with the target rendering system. The version used in this thesis is 3.7.0. Vega oper¬ 
ates on Pentium III platforms or superior, and requires 64 MB of RAM (minimum). A 
graphics accelerator card with at least 32 MB is required to increase performance. It runs 
in Windows 2000, Windows NT and SGI Irix. 

The goal of Vega is to provide a system to the developer that can be used to con¬ 
struct a simulation application quickly, with the highest performance possible. To ac¬ 
complish this, Vega includes Eynx, a graphical interface that allows the user to construct 
the bulk of the application in a point and click environment. Eynx generates an Applica¬ 
tion Definition Eile (ADE) that is used to configure the application. Eynx also allows the 
user to preview the application and make adjustments without leaving the Eynx environ¬ 
ment. An ADE may be used to configure a stand-alone executable, or in conjunction with 
user-supplied code, in order to create complex visual applications. 

2, The Vega Marine Module 

The Vega Marine module provides the capability for Vega to simulate moving 
ocean waves along with marine special effects. The ocean is simulated as a dynamic 
component extending to a radius about an observer, and a static component extending 
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further. The dynamie oeean ereates a polygonal approximation to a wave field defined by 
the superposition of 10 sinusoidal components. All components propagate in the same 
(user-specified) direction. A default set of sinusoidal components is included in the Vega 
Marine module along with a set of sea states for which the default components are de¬ 
rived using energy spectrum analysis [5]. 

The virtual waves are created by the Vega Marine module using a sum of 10 sinu¬ 
soids and a user defined set of non-harmonic frequencies selected to produce a realistic- 
appearing wave field over a range of sea states. The instantaneous wave elevation is a 
space-time function computed as follows: 

10 

4'(x,y,0 = f + 0.112//^^ cos(k, {xsinz + y cos z) + + (p^) (4.1) 

f=l 

where 

T = tide, the average height above mean sea level (meters), 

//s = fundamental wave height for the current sea state (meters), 
k\ = (f2i cOs)^ / g, deep-water dispersion relation of the component, 
g = gravity acceleration (9.81 m/s ), 

X = wave direction (0 to In), 

= non-dimensional angular frequency of the component, 

©s = 27 i/J's, angular frequency of sea state, 

Ts = the modal period of the sea state, 
t = time (s), and 

(pi = phase angle of the /* component. 

Vega Marine provides application program interface (API) allowing the host ap¬ 
plication to get and set the parameters associated with the virtual representation of the 
dynamic ocean. This, in conjunction with a user-furnished dynamics model, allows the 
host to produce the correlated motion of the ship 


14 



3, Special Effects 

Vega Marine provides many of the marine speeial effeets required by today’s so- 
phistieated maritime training applieations. It ineludes astern wakes, bow waves, eddies, 
buoys, eonstant-length lines, eonstant-tension lines, erew overboard, hoisted flags, knuek- 
les, whiteeaps, flotsam, foam, log and surf The models of these speeial effeets are not 
open to users, but it is possible to adjust some internal parameters to aehieve more realis- 
tie results. This simulator uses astern wake, bow wave, buoys and eddies. 

Furthermore, Vega provides all the resourees to eontrol the lightning of the seene, 
ineluding fog effeets. An additional resouree that is used to inerease the realism of the 
simulation is the ship’s smoke. The smoke follows the wind relative direetion, eonsider- 
ing the speed of the ship. 

4, Vega Coordinate System 

Vega provides three types of eoordinate systems; flat, spherieal and elliptieal. 
This work sets Vega in the flat eoordinate system, whieh means that the three axes Xy, Yy 
and Zy, are oriented as east, north, and up respeetively, as depieted in Figure 2. 

Zv 



Sinee the eoordinate systems used by models may vary, some homogeneous trans¬ 
formations are required to adjust the output parameters to the eorreet rendering frame. 
Following the deseription of eaeh model used, the respeetive required homogeneous 
transformation is deseribed. 
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This chapter described the main features of the Vega system that is the basis of 
building the Virtual Environment module. Following the description of the system mod¬ 
ules, the next chapter describes the Ship Motion module in detail. Initially, a theoretical 
background is introduced, and then all models related to the motion of the ship are com¬ 
pletely explained. 
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V. THE SHIP MOTION MODULE 


This chapter first presents some important coneepts regarding a ship’s hydrody- 
namies in order make it easier to understand the simulation. The ehapter then describes 
the models used in this simulation, starting with the Maneuvering Model. Two wind¬ 
generated wave models are used. The Wave Model I is a six degrees-of-freedom (DOF) 
wave model that works deeoupled from the Maneuvering Model. The Wave Model II is 
eoupled to the Maneuvering Model and adds a random component to the roll moment. 
The description of the Wind and Ocean Current Models completes this chapter. 

A, THEORETICAL BACKGROUND 

1. Degrees of Freedom and Notations 

The motion of marine vehieles in six DOF is eonsidered since six independent 
eoordinates are neeessary to determine the position and orientation of a rigid body. The 
deseription of each DOF and the respective notation are shown in Table 2. 


DOF 

index 

Description 

Forces and 
moments 

Linear and an¬ 
gular veloc. 

Positions and 
Euler angles 

1 

Motions in the x-direetion (surge) 

X 

U 

X 

2 

Motions in the y-direetion (sway) 

Y 

V 

y 

3 

Motions in the z-direetion (heave) 

Z 

w 

z 

4 

Rotations in the x-direetion (roll) 

K 

p 

(p 

5 

Rotations in the y-direetion (piteh) 

M 

q 

e 

6 

Rotations in the z-direction (yaw) 

N 

r 

¥ 


Table 2. DOF Description and Notation. 


The first three coordinates and their time derivatives correspond to the position 
and translational motion along the x-, y- and z-axis, while the last three coordinates and 
time derivatives are used to describe orientation and rotational motion. 

Using the above notation, the general motion of a marine vehicle in six DOF can 
be described by the following vectors (bold variables are vector variables): 
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(5.1) 


^ = ^i=[x,y,zf, x\^ =[(!), e, y/f, 

D = [D[','i)2f, 'l)i=[w,v,wf, v^=[p,q,rf, 

T = {T\,xlf, T,=[X,7,Zf, x,={K,M,m". 


2, Coordinate Frames 

When analyzing the motion of marine vehicles in six DOF, it is fitting to define 
two coordinate frames, as shown in Figure 3, from Ref. 3. 



X 


o 


X 


Figure 3. Maneuvering Model Coordinate System. 

The moving coordinate frame Xo Yq Zq is fixed to the vehicle and is called the 
body-fixed reference frame. The origin O of this frame is usually chosen to coincide with 
the center of gravity (CG) of the ship. For marine vehicles, the body axes Xq, Yq and Zq 
coincide with the principal axes of inertia. The motion of the body-fixed frame is de¬ 
scribed relative to an inertial reference frame XY Z. 

3, Ship Motion Equations 

Representing the motion equations in the Cartesian system of coordinates (body- 
fixed reference frame) and defining xg, yc and zg as the position of the ship’s CG, the 
well known motion equations of a rigid body [3] are giving by the following; 
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(5.2) 


X = m\u + qw-rv + XQ{q^ +r^) +yQ{pq-r) +ZQ{rp + q)'\ 

Y = m[v + ru-pw-y^ir^ + p^) + Z(.{qr-p) + X(;{qp + r)] (5.3) 

Z = m[w + pv-qu-z^{p^ +q^) + x^ {rp -q) + yQ {rq + p)] (5.4) 

K = Ij^p + (I^ -ly)qr + m[yQ(w + pv-qu)-ZQ(u + m-pw)\ (5.5) 

M =Iyq + {I„ -1 2 )rp + m[zQ{u + qw - rv) - Xq{w + pv-qu)\ (5.6) 

N = I^r + {Iy -I^)pq + m[Xg{v + m-pw)-yQ{u + qw-rv)] (5.7) 


where m is the ship’s mass; u,v,w,u,v,w represent the linear veloeities and aeeelerations 
in the xo, yo and zo direetions; and p, q, r, p, q, r represent the angular veloeities and ae- 
celerations related to the axes xq, yo and zq. Ix, h and Iz are the inertia tensors related to 
same axes of the body-fixed frame. The forees and moments X, Y, Z, K, M and N repre¬ 
sent the resultants of all external aetions over the ship’s body. 

These equations can be expressed in a more compact form, i.e., the vectorial rep¬ 
resentation of the 6-DOF rigid-body equations of motion equations: 

^RB^ ^ ^~ RB (5-8) 

where Mrb is the inertial matrix, Crb is the matrix of Coriolis and centripetal terms and 
Trb is the vector of forces and moments. 

In basic hydrodynamics, it is common to assume that the hydrodynamic forces 
and moments on a rigid body can be linearly superimposed by considering the radiation- 
induced forces and diffraction forces [11]. The radiation-induced forces include added 
mass due to the inertia of the surrounding fluid, the radiation-induced potential damping 
due to the energy carried away by generated surface waves, and restoring forces due to 
Archimedes Law (weight and buoyancy). 

Therefore, Eq. 5.8 can be augmented, yielding a more complete form to the vecto¬ 
rial representation of the 6-DOF motion equations: 
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Ml) + C(d)d + Z)(d)d + g(ti) = T 


(5.9) 


Ti = j(n)v (5.10) 

where M is the inertial matrix (including added mass), C(v) is the matrix of Coriolis and 
centripetal terms (including added mass), D(v) is the damping matrix, g(x\) is the vector 
of gravitational forces and moments and t is the vector of the propulsion forces and mo¬ 
ments. 


The concept of added mass is usually misunderstood to be a finite amount of wa¬ 
ter connected to the vehicle such that the vehicle and the fluid represent a new system 
with mass larger than the original system. Added (virtual) mass should be understood as 
pressure-induced forces and moments due to a forced harmonic motion of the body, 
which is proportional to the acceleration of the body [3]. 

In Eq. 5.10, J is the transformation matrix used to represent the position and ori¬ 
entation vector in the Earth-fixed frame. Equation 5.10 is rewritten below with all de¬ 
tails of J; 


X 


y 


z 




e 


y. 



cy/ 

sy/cO 

-sO 


-sy/ccj) + cy/s6s(j) 
cy/ccj) + s(j)s6sy/ 
c6s(j) 


0 


3x3 


sy/s(j) + cy/c(j)s6 
-cy/scj) + s6sy/c(j) 


^3x3 


C0C(j) 

1 

S(l)tO 

c(j)tO 


0 

C<j) 

-S(j) 


0 

scj)! cO 

C(j)l C0 


where s = sin, c = cos and t = tan. 


4. Non-Dimensional Parameters 

When using hydrodynamic models, it is often suitable to normalize the ship steer¬ 
ing equations of motion such that the model parameters can be treated as constants with 
respect to the instantaneous speed U. There are many normalization forms, but the most 
commonly used for this kind of application is the Prime-system of SNAME [14]. Table 3 
shows the Prime-system factors to be applied depending on the unit of the variables of in¬ 
terest. 
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Unit 

Factor 

Unit 

Factor 

Length 

L 

Linear velocity 

U 

Mass 

pL^/2 

Angular velocity 

U/L 

Inertia moment 

pL^/2 

Linear acceleration 

U^/L 

Time 

L/U 

Angular acceleration 

U^/L^ 

Area 


Force 

pU^L^/2 

Position 

L 

Moment 

pU^L^/2 


Table 3. Non-Dimensional Prime-System Faetors. 


This thesis uses the ship’s instantaneous speed U, the length L = Lpp (between the 
fore and aft perpendiculars) and the salt-water density p as normalization variables. The 
notation used to indicate non-dimensional variables is the “prime” symbol ('). For exam¬ 
ple, v' = v/f/ is the normalized (non-dimensional) form of the velocity v. 


5. Hydrodynamics Derivatives 

An important step in the development of the maneuvering models is to expand the 
forces and moments in Taylor’s Series. In this way, the nonlinear terms show up as pow¬ 
ers of the independent variables, establishing a polynomial equation. To use the expan¬ 
sion, the function and its derivatives need to be continuous and finite in the region of val¬ 
ues of the variables under consideration. The precision of the model depends upon the 
point where the series is truncated. Most models expand the series up to the third power. 
Some terms can be neglected due to its low significance or because they are null due to 
the symmetric geometry of the ship. The equation below shows an example of the sway 
force Y in this form of representation; 

Y = Yp> + Y/ + Y^p + Y^(j) + Y^y + Yy + Y^r + Y^y ^ 

++Y,,/(!> + 


where the terms that characterize the hydrodynamic derivatives are represented as 

, and so forth, evaluated at equilib- 


T. = 


fQY'^ 


d(j) 


Y = 

9 rrr 


^ dy 

y^Jr 


Y = 

9 vvr 


f Q^Y ^ 


d vdr 
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rium conditions. The initial condition of motion equilibrium is chosen as straight ahead 
motion at a constant speed [14], These terms are also known as hydrodynamic coeffi¬ 
cients. 

A large number of experimental methods are available to determine forces and 
moments associated with variations in linear and angular velocity and acceleration. Typi¬ 
cal facilities are the rotation arm, the free oscillator, the forced oscillator, the curved-flow 
tunnel, the curved models in a straight flow facility and the Planar Motion Mechanism 
(PMM) technique. Nevertheless, it is difficult to determine all the hydrodynamic coeffi¬ 
cients for an ocean vehicle. It is necessary to know these coefficients with reasonable ac¬ 
curacy to obtain a good model of the vehicle [3]. Thus, it is important to notice that the 
quality of the models relies on experimental databases. 

6. The Environment Disturbances 

Expanding the nonlinear dynamic equation of motion (Eq. 5.9), it becomes: 

+ (u) + (u)]u + [Dp (u) + Dp (u)]u + g(ti) = t (5.13) 

where Mrb and Ma are the mass and added mass matrices, Crb and Ca are the Coriolis 
and centripetal terms matrices for the rigid-body and added mass, and Dr and Dy are the 
radiation-induced potential damping and viscous damping matrices, respectively. 

Applying linearization about zero angular velocity and u = uo, v = vq and w = wo, 
the Coriolis and centripetal terms collapse to zero and Eq. 5.13 can be written as follows: 

[Mpp+MJv + [Np + Np]v + Gr\ = T (5.14) 

where N, = and G = are three eonstant 

dr] 

matrices. The principle of superposition, one of the properties of linear systems, suggests 
that the environment disturbances can be added to the right hand side of Eq. 5.14 to yield: 

Wrb +MJv + [Np + Np]v + G^ = ^ + . (5.15) 

Eet the ocean current velocity vector be defined by 0,0,0]^, where 

the three last fluid velocity components are zero, assuming irrotational fluid. Assuming 
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also that the vehicle is neutrally buoyant, the mass is homogenously distributed and 
=i)-i)^can be interpreted as the relative velocity vector, the linear equations of mo¬ 
tion can be combined to yield the following: 

RB~^ ^ ]'*^r + ^ ^ ^ wave ^ wind * (5.16) 

An extension to the nonlinear case proposed by Fossen [16] considers that the fol¬ 
lowing approximation (Eq. 5.17) can be used by ships in the case of treating slowly vary¬ 
ing currents in terms of the relative velocity: 

[Mub + ^ (u J (u Jjt), + Z)(u Ju, + g(tl) = t + + T . (5.17) 

B, THE MANEUVERING MODEL 

An important decision in this study is concerned with the type of ships that would 
be used. The choice was based on the availability of pre-developed dynamic models and 
the associated environmental disturbance models. It is preferable that the ship belongs to 
a very generic class of merchant ships, so that most of the needed information is not clas¬ 
sified. 

Searching in the literature, a very suitable ship model was found that was written 
in the Matlab language. It is a 4-DOF nonlinear rolling coupled steering model for a high¬ 
speed container ship [3]. In general, for such surface ships, it is common to find models 
that neglect the heave, pitch and roll modes under the assumption that these motion vari¬ 
ables are small. Fortunately, this model includes the roll motion variable that renders the 
final results more realistic. 

The principal limitation of this model is that it only allows forward engine orders. 
It not only constrains backward maneuvers but also impedes some crash stop maneuvers. 

1. Container Ship Main Features 

The container ship is described by the set of data in Table 4. 
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Description 

Symbol 

Value 

Length 

L 

175.00 m 

Breadth 

B 

25.40 m 

Draft 

fore 

dp 

8.00 m 

aft 

dA 

9.00 m 

mean 

d 

8.50 m 

Displacement volume 

n 

21,222.00 m^ 

Height from keel to transverse metacenter 

KM 

10.39 m 

Height from keel to center of buoyancy 

KB 

4.6154 m 

Block coefficient 

Cb 

0.559 

Rudder Area 

Ar 

33.0376 m^ 

Aspect Ratio 

A 

1.8219 

Propeller Diameter 

D 

6.533 m 


Table 4. Container Ship Main Features. 


2, Nonlinear Equations of Motion 

A mathematieal model for a single-serew high-speed container ship in surge, 
sway, roll and yaw was presented by Son and Nomoto [12], In this model, the rigid-body 
dynamics including the contribution from the hydrodynamic added mass derivatives are 
presented as follows: 


X' = (m' + m'^)u -(m' + m'y)v'r' 

(5.18) 

Y' = (m' + m'y)v' + (m' + m'Ju'r' + m'^a'^r' - m'yr^,p' 

(5.19) 

K' = ii:+ J'x)p' - m'// - m'/y/ + W'GM'</>' 

(5.20) 

N' = (/;+ j'y + m'yy + xy. 

(5.21) 


All the equations are in the non-dimensional form, denoted by the symbol ('). The 
variables m'x, m'y, J'z and Jx denote the added mass and added moment of inertia in the x- 
and y-directions about z- and x-axes, respectively. The center of added mass for m'y is de¬ 
noted by a'y (x-coordinate) while I'x and I’y are the added mass z coordinates of m'x and 
m'y, respectively. Notice that these equations derive from the standard motion equations 
of a rigid body (Eq. 5.2-7) with heave and pitch neglected. 
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The metacentric restoring moment in roll was added in the K' expression (Eq. 
5.20), where W is the weight of the water displaced by the ship hull and GM' is the 
transverse metacentric height. 


The hydrodynamic forces and moments are defined in equations below: 

X'= x'{u ')+(1 - t)r{j) + xy/ + xy^ + xy^+xyy c^f', sin (5.22) 


7' = ry + Yy + ry + ry ++ o'V'+CvV'^ 

+ O'V' + + K/y + Y'/y + (1 + a, )F; cos 


K' = K'y+K'y+K'y+ry++ xy^ + xyy + xyy 
+ xyy + xyy+xyy+xyy - (i+«, )z;f; cos 


N'=N'y+N'/+N'y + N'y + xy^ + xy^ + xyy+xyy 

+V'+ yy<t >^+V'+ K6/<t >^+ (4 +cos 5' 


(5.23) 


(5.24) 


(5.25) 


-y 

where X(m) is a velocity-dependent damping function, Q.g.,X(u) = Xuuu' . 


The hull parameters used in this model are presented and identified in Table 5. All 
parameters are in their non-dimensional form. 


Parameter 

Description 

Value 

m' 

Mass of the ship 

0.00792 

m’x 

Added mass in the x direction 

0.000238 

m'y 

Added mass in the y direction 

0.007049 

I’x 

Moment of inertia about x axis 

0.0000176 

Jx 

Added moment of inertia about x axis 

0.0000034 


Moment of inertia about z axis 

0.0000176 

J. 

Added moment of inertia about z axis 

0.0000034 

a'y 

Center of added mass m'y (x-coordinate) 

0.05 

I'x 

Added mass z-coordinate of m'x 

0.0313 

I'y 

Added mass z-coordinate of m'y 

0.0313 

X’uu 

Hydrodynamic derivative d^X!du^ 

-0.0004226 

X'„ 

Hydrodynamic derivative d^X / dvdr 

-0.00311 

X'xx 

Hydrodynamic derivative d'^X / dv^ 

-0.00386 
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Parameter 

Description 

Value 

Xrr 

Hydrodynamic derivative X / dr^ 

0.00020 

^ (pip 

Hydrodynamic derivative d^X / dcj)^ 

-0.00020 

Y, 

Hydrodynamie derivative dY / dv 

-0.0116 

Yr 

Hydrodynamie derivative dY / dr 

0.00242 

Yp 

Hydrodynamie derivative dY / dp 

0.0 

Y, 

Hydrodynamie derivative dY / dcj) 

-0.000063 

V 

1 vvv 

Hydrodynamie derivative d^Y / dv^ 

-0.109 

V 

± rrr 

Hydrodynamie derivative d^Y / dr^ 

0.00177 

V 

1 rvv 

Hydrodynamie derivative d^Y / drdv^ 

0.0214 

V 

1 rrv 

Hydrodynamic derivative d^Y / dr^dv 

-0.0405 

V 

vv^ 

Hydrodynamie derivative d^Y / dv^dcj) 

0.04605 

Y 

^ Vf(p 

Hydrodynamie derivative d^Y / dvdcj)^ 

0.00304 

V 

^ rr(p 

Hydrodynamie derivative d^Y / dr^dcj) 

0.009325 

V 

J- r(p(p 

Hydrodynamie derivative d^Y / drdcj)^ 

-0.001368 

iV'v 

Hydrodynamie derivative dN / dv 

-0.0038545 

N'r 

Hydrodynamie derivative dN / dr 

-0.00222 

N’p 

Hydrodynamie derivative dN / dp 

0.000213 

N'^ 

Hydrodynamie derivative dN / dcj) 

-0.0001424 

iV'vw 

Hydrodynamie derivative d^N / dv^ 

0.001492 

N’rrr 

Hydrodynamie derivative d^N / dr^ 

-0.00229 

N’rw 

Hydrodynamie derivative d^N / drdv^ 

-0.0424 

N’rrv 

Hydrodynamic derivative d^N / dr^dv 

0.00156 

N' wip 

Hydrodynamie derivative d^N / dv^dcj) 

-0.019058 

v(p(p 

Hydrodynamie derivative d^N / dvdcj)^ 

-0.0053766 

N’ rrf 

Hydrodynamie derivative d^N / dr^dcj) 

-0.0038592 

r(p<p 

Hydrodynamie derivative d^N / drdcj)^ 

0.0024195 

K'p 

Hydrodynamic derivative dK / dp 

-0.0000075 

K'^ 

Hydrodynamie derivative dK / d(j) 

-0.000021 

iV vw 

Hydrodynamie derivative d^K / 

0.002843 

iV rrr 

Hydrodynamie derivative d^K / dr^ 

-0.0000462 

iV rvy 

Hydrodynamie derivative d^K / drdxY 

-0.000558 

iV rrv 

Hydrodynamie derivative d^K / 

0.0010565 

vv^ 

Hydrodynamie derivative d^K / dv^dcj) 

-0.0012012 

V(p(p 

Hydrodynamic derivative d^K / dvdcj)^ 

-0.0000793 
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Parameter 

Description 

Value 

iV rr(p 

Hydrodynamie derivative d^K / dr^dcj) 

-0.000243 

rqxp 

Hydrodynamic derivative d^K / drdcj)^ 

0.00003569 


Table 5. Hull Parameters in the Maneuvering Model. 


The rudder normal foree is resolved as follows: 


F' = 

^ N 


6.13A 4^ 2 
A + 2.25 "" 


+ v^ lsm«^ 


(5.26) 


where u'r and v'r are the surge and sway eomponents of the ineident flow veloeity in the 
rudder, and ur is the ineidenee flow angle, defined by Eqs. 5.27, 5.28 and 5.29, respee- 
tively: 


U„ = UpS. iH-^ 


v; = yv + c „/+ 


a,^=5 + tan 




J 


(5.27) 

(5.28) 

(5.29) 


Equations 5.30 and 5.31 define/, the advanee number, and u'p, respeetively: 


_ UpU 

nD 


(5.30) 


u’p = eos v'[(l - %) + r{(v' + x’/f + c^y + c^/}]. (5.31) 

The Bilinear Thruster Model is defined by Eq. 5.32, where T is the thrust devel¬ 
oped by the propeller and Kr (J) is the thrust eoeffieient: 


T = lpD^Kp{J)\n\n 

(5.32) 

Ar(/) = 0.527-0.455/. 

(5.33) 


The propeller and rudder parameters used in this model are presented and identi¬ 
fied in Table 6. 
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Parameter 

Description 

Value 

t 

thrust deduction number 

0.175 

Wp 

wake fraction number 

0.184 

x’r 

x-coordinate of point on which the rudder lateral force Yg acts 

-0.5 

z'r 

z-coordinate of point on which the rudder lateral force acts 

-0.5 

x'p 

x-coordinate of propeller position 

-0.526 

Qh 

rudder to hull interaction coefficient 

0.237 

x’h 

x-coordinate of point on which normal force Fm acts 

-0.48 

Crx 

constant 

0.71 

T 

constant 

1.09 

s 

constant 

0.921 

k 

constant 

0.631 

y 

flow rectification coefficient 

0.088 

Cpv 

propeller flow rectification coefficient 

0.0 

^pr 

propeller flow rectification coefficient 

0.0 

^Rr 

rudder wake coefficient 

-0.156 

Cprrr 

Rudder wake coefficient 

-0.275 

^Rrrv 

Rudder wake coefficient 

1.96 


Table 6. Propeller and Rudder Parameters in the Maneuvering Model. 


3, The Model Representation 

The maneuvering model is basically represented by Eqs. 5.9 and 5.10. Some sim¬ 
plifications were carried out and adjustments were made to fit only 4-DOF. As a result, 
the model’s equations are; 

Mv = t (5.34) 

ti = J(tl)n (5.35) 


The matrix of mass, added mass and inertia M is derived from Eqs. 5.18-21 and 
expressed as 


m' + m[ 0 0 

0 m' + m' -m'T, 

0 

0 0 

y } 


myUy 

0 

H+Jl 


(5.36) 
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Since M is non-singular, the first set of state derivatives can be written as: 

= M-' . (5.37) 

A 

_N'_ 

In order to simplify the notation and taking advantage of its symmetry, the matrix 
M is re-written as: 

m, j 0 0 0 

0 ^22 ^32 W 42 

0 ^32 ^33 0 

0 ^42 0 ^44 

After computing the inverse of matrix M, Eq. 5.37 is explicitly written as follows: 

(l/m,i)X' 

(-m32m447' + m22«44X'-m42X' + m32m42X')/(m22^33^44 “^32^44 “^ 42 ^ 33 ) 

+ m22in^2^' + m22tn22N' - ^ 3 ^A')/(m22m33m44 - «3\«44 - ml2fn22) 

Returning to the dimensional form, the factor jh is applied to the linear veloc¬ 
ity derivatives and the faetor U^jl} is applied to the rotational veloeity derivatives, 
yielding the final form: 

m] \u^Il 0 0 0 1 fii' 

V _ 0 U^jL 0 0 v' 

p ~ Q 0 U^jL^ 0 p 

f J L 0 0 0 U^Il^\ Ir 

From the 4-DOF adjusted Eq 5.35, the other set of state derivatives ean be written 
as below: 


(5.40) 




(5.38) 
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(5.41) 


• t 

X 


cost// 

-sin^cos^zi 
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sini// 

cosy/cos (/) 

0 
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t 

V 



0 

0 

1 
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t 

p 



0 

0 

0 

COS(/) 


r' 


Again, to return to the dimensional form, the faetor U is applied to the position 
derivatives and the faetor UjL is applied to the angle derivatives, as shown below: 
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■[/ 0 0 0 “ 
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¥' 


(5.42) 


4, Rudder Saturation and Dynamics 

The rudder dynamics is based on a first order model suggested by Van Ameron- 
gen to represent the steering machine [15]. This steering machine model is depicted in 
Figure 4. Generally, the rudder angle and rudder rate limits are typically in the ranges of 
-35“ < < 35“ and 2.33“/s < ^nax < ’ respectively. In this model, according to the 


characteristics of the container ship, the rudder angle and rudder rate limits were set at 
20“ deg and 3“/s, respectively. The parameter Sc is the commanded rudder angle. 



Figure 4. The Steering Machine Model. 

5, Shaft Speed Saturation and Dynamics 

The shaft dynamics used in this model is also a first order system. The behavior of 
this system depends upon the current shaft speed value. There are two ranges of opera¬ 
tion, which determine the value of the parameter Tm to be used. The first one is for speed 
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values lesser or equal to 20 rpm, with being a eonstant value equal to 18.83. The see- 
ond is for values greater than 20 rpm, with Tm being a funetion of n, defined by 
= 5.65/n . Figure 5 shows the shaft model. 
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Figure 5. The Shaft Model. 

6, Interaction with Other Modules 

The Maneuvering Model module has basieally two input parameters, the eom- 
manded shaft speed {ric) and the eommanded rudder angle (<5^). These parameters eome 
from the User’s I/O sub-module of the Virtual Environment module. Indeed, the shaft 
speed order is obtained from the engine order command. Each command corresponds to a 
shaft speed; Stop = 0 rpm; /s Ahead = 53 rpm; /s Ahead =107 rpm; and Eull Ahead = 
160 rpm. The rudder angle order range accepted by this model is from -20° to 20°. 

The Maneuvering Model module provides ten state derivative variables: surge and 
sway linear accelerations, roll and yaw angular accelerations, surge and sway linear ve¬ 
locities, roll and yaw angular velocities, the rudder angular velocity and the shaft angular 
acceleration. In order to obtain the parameters of interest, it is necessary to integrate these 
state derivatives over the elapsed time between the renderings of two frames by the 
Graphics module. This amount of time. At, is measured at each frame-rendering iteration 
by the Vega’s function vgGetDeltaFrameTimeQ. The value is passed to the thread that 
treats the models to perform the integration. There are many methods to numerically in¬ 
tegrate functions [13]. The method used in this thesis was the Euler’s integration method 
due to its simplicity. It worked well in this application because the sampling rate is high 
(about 60 Hz) compared to the natural frequency of the model. The state derivatives are 
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multiplied by At and added to their previous values to generate the state variables that are 
passed to the Virtual Environment module to define the new position and orientation of 
the ship for the next frame rendering, as below: 

(5.43) 

The aetual values of the rudder angle and the shaft speed are passed to the Virtual 
Environment module as well in order to be presented on the terminal sereen. 

The last step before passing the new position and orientation to the Virtual Envi¬ 
ronment module is to apply a eoordinate transformation to eonvert the results from the 
Maneuvering Model Coordinate System to Vega Coordinate System. Eigure 6 shows the 
two frames involved in this operation. 



Eigure 6. Maneuvering Model to Vega Coordinate System Conversion. 


The eonversion is aehieved by applying a transformation matrix to every position 
generated by the Maneuvering Model. It must rotate 180 degrees by the Xm axis and then 
rotate -90 degrees by the Zm axis. Equation 5.44 shows the homogeneous transformation 
required. 
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(5.44) 


The parameters related to rotational motions (roll and yaw) follow the right hand rule in 
both frames. 
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The Wave Model II and the Wind Model are coupled to the Maneuvering Model. 
The interactions among them are discussed in their respective sections. 


7, Validation of the Maneuvering Model 

The idea of this validation section is not to validate the original Matlab code since 
the model was obtained from a very reliable source. The part that required validation is 
the C++ implementation of the same model. However, the conversion from Matlab to 
C++ was very successful. No numerical issues were observed, so that the results obtained 
in a comparative test were very much the same. Figure 7 shows the surge and sway pa¬ 
rameters from the two sources. It is seen that the two output plots (from Matlab original 
code and from C++ implementation) are perfectly the same, one overlapping the other. 


Comparing the surge velocity 




Figure 7. Matlab Code vs. C++ Implementation. 


Instead of continuing the presentation of these coincident plots, the rest of this 
validation section was used to evaluate some performance features of the maneuvering 
model. At first, from the acceleration/deceleration test, it was verified that the time spent 
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to go from zero speed to the maximum speed (engine order = Full Ahead) was 355 sec¬ 
onds, approximately 6 minutes. The time to do the inverse procedure (engine order = 
Stop) was greater than 30 minutes. Figure 8a shows the acceleration/deceleration plot. 


(a) Acceleration / Deceleration 




Figure 8b demonstrates that the shaft spends approximately 20 s to reach its 
maximum speed and approximately 85 s to completely stop. 

From the rudder dynamics results in Figure 9, it is seen that the time spent to go 
from the neutral position (0°) to the maximum angle allowed by the model (± 20°) was 
approximately 12 s. The time to the maximum excursion, from -20° to 20° was approxi¬ 
mately 20 s. Since the rudder angle rate did not exceed the limit of 3°/s, i.e., the normal 
rate for that kind of ship, the result is as expected. 

Figure 10 shows the behavior of the surge and sway velocities during maneuver¬ 
ing. As expected, the surge velocity decreases and the sway velocity increases in the op- 
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posite direction of the movement during the turnings, according to the inertial reaction. It 
is also seen from this plot that the surge and sway velocities have a tendency to return to 
their initial values after the end of all maneuvering. 


Rudder Dynamics 



Figure 9. Rudder Dynamics. 


Surge and Sway Velocities 



Surge and Sway Velocities. 
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Figure 10. 











































Other feature observed in this evaluation process was the roll and yaw behaviors 
during the maneuvers. From Figure 11, it is seen that roll angle has the expected behav¬ 
ior, i.e., the tilt of the ship is in the opposite direction of the turning, according to the in¬ 
ertia. The yaw angle also had the expected performance. The plot shows that, from six to 
nine minutes, the yaw angle varied approximately 200°, implying a turning rate of 
0.6677s. That rate is perfectly acceptable for that class of ships. 
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Yaw Angle 
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Figure 11. Roll and Yaw Angles. 

Two more classic tests were done; the turning circle and the zig-zag maneuvers. 
The turning circle test was initiated with the ship at a speed of 21.6 knots. At an arbitrary 
point, a 20° rudder angle turn to starboard started. The parameters evaluated in this test 
were the steady turning radius (0.245 nm) and the time spent in complete maneuver (ap¬ 
proximately 500 seconds). Compared with data from similar ships, those values were 
considered perfectly valid for this kind of application. Figure 12 shows the turning circle 
test plot. 



Roll Angle 
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Turning Circle 



Figure 12. Turning Circle. 

The last evaluation test was the Kempf s zig-zag test. This test was first proposed 
by the German scientist Gunther Kempf in 1932 [19]. The zig-zag time response is ob¬ 
tained by commanding the rudder 20° starboard from an initially straight course. The 
rudder setting is kept constant until the heading is changed 20°, then the rudder is re¬ 
versed 20° to port. Again, this rudder setting is maintained until the ship heading has 
reached 20° in the opposite direction. The test continues until a total of five rudder 
changes have been completed. 

The zig-zag test is usually referred to as 20°-20° maneuver, where the first angle 
refers to the actual rudder settings and the second angle denotes how much the heading 
angle should change before the rudder is reversed. For large ships, it is recommended to 
use of a 10°-10° or a 20°-10° maneuver to reduce the time and the space required. In this 
thesis, the original 20°-20° maneuver was used in order to be possible to compare the re¬ 
sults with the available data. 
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The test started at the eonstant speed of approximately 11 knots. The elapsed time 
to reaeh the first overshoot was 82 s. The average overshoot angle was approximately 
6.5°, and the elapsed time to eomplete the first eyele was 295 s. The final speed was ap¬ 
proximately 8 knots. Figure 13 shows the end result of the test. 

The results were very eompatible with what eould be expeeted from a large eon- 
tainer ship behavior. The visual validation corroborated these results. The final evaluation 
of the maneuvering model was completely satisfactory. 


Kempf's Zig-Zag Maneuver 



Figure 13. Kempfs Zig-Zag Maneuver. 

C. THE WAVE MODEL I 

The Wave Model 1 is a 6-DOF motion model, ft is decoupled from the Maneuver¬ 
ing Model so that its results do not influence ship dynamics. As mentioned previously, in 
order to achieve good results in terms of realism with a low computational cost, an inter¬ 
polation technique was implemented in the wave model simulation. This technique al¬ 
lows the use of the results of a complex wave model in a very computationally efficient 
way, by interpolating the desired parameters from tables previously generated. In order to 
make this possible, some assumptions were made; 1) the ship’s speed varies from 0 to 60 
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ft/s; 2) the wave heading (WH) varies from 0° to 360°, as shown in Figure 14; 3) the 
wavelength of the fundamental eomponent of the oeean model varies from 10 to 1460 ft; 
and 4) the ship has a eonstant displacement during the entire simulation. Each table pro¬ 
vides, for a given (speed, wave heading angle) pair, the values of the output parameters 
related to the set of pre-defined wavelength values of the fundamental component of the 
current sea state. The outputs are surge, heave, pitch, sway, roll and yaw. Each of them 
has amplitude and phase parameters. The computer program described below generates 
these tables. 


Beam Sea 
WH = 90“ 



Eigure 14. Wave Heading Notation. 

1. General Description 

The computer program called SHIPMO.BM, by Robert E. Beck and Armin W. 
Troesch [6], was developed to predict ship motions in six degrees of freedom and five 
components of the shear and bending moment distribution. It is based on the strip theory 
of Salvesen, Tuck and Ealtinsen [7] and it was extended by Beck [8] to include the surge 
degree of freedom. The program computes both regular and irregular wave results. 

This program, originally written in EORTRAN, requires two input data files. The 
first, SHIPMO.IN, contains basically three types of information: some internal parameter 
settings to define the desired operation; information about the ship’s structure and the de- 
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sired output ranges, in terms of ship’s speed, wave heading angle and wavelength. The 
second, G2S1.3, contains terms that are needed in shallow water but not for deep-water 
calculations, assumed for this thesis, but must be present for the program to run. 

The program generates output files in different formats. The file ADMASS con¬ 
tains hydrodynamic added mass and damping results. The fde matdata is an output fde 
with data in Matlab-ready format, while POT and RAO are files that contain results that 
could be used in subsequent runs of the program so that it runs faster. The fde VISC con¬ 
tains viscous roll damping results, and SHIPMO.OUT is the main text-based results file. 

The SHIPMO.OUT file provides non-dimensional values in both vertical and 
horizontal plane responses. These values are related to surge, heave and pitch parameters 
in the vertical plane and sway, roll and yaw parameters in the horizontal plane. Each pa¬ 
rameter has amplitude (non-dimensional) and phase (degrees). This phase angle repre¬ 
sents the difference of phase between each parameter and the encounter wave created by 
the ship motion based on the Doppler effect. In order to convert the outputs to the dimen¬ 
sional form, the linear motion parameters (surge, sway and heave) are multiplied by the 
current wave amplitude in the ship’s CG. The rotational outputs (pitch, roll and yaw) are 
multiplied by the same wave amplitude times the wave number. The wave number is de¬ 
fined as2;r//l, where 2 is the wavelength of the fundamental component of the current 
sea state. 

2. Wave Model Coordinate System 

The coordinate system adopted by this model is depicted in Figure 15. It has the 
X-axis as the longitudinal axis of the ship (directed from aft to fore), the T-axis as the 
transverse axis (directed to port) and the Z-axis as the vertical axis (directed up). 
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3, Wave Model Tables 

In order to obtain the set of tables to exeeute the interpolation proeess in the simu¬ 
lation program, two distinet steps are required; 1) run the program to gen¬ 

erate the output file SHIPMO.OUT with the tables, and 2) convert this output file to a 
C++ type waveModelTables.h file. This file contains the tables in two-dimensional array 
format turning the tables readable to the interpolation program. 

a. Obtaining SHIPMO. OUT File 

In order to perform the first task, the wave model’s input file SHIPMO.IN 
must be prepared. The internal parameters were set as shown in Table 7. 


Settings 

Description 

1A = 0 

No even spacing 

1B = 0 

Zero weight curve ordinates 

1C = 0 

Units on displacement = metric or long tons 

1D = 0 

Regular waves model 

1E= 1 

Degrees of freedom = both vertical and horizontal plane 

1F = 0 

Eong crested seas 

1G= 1 

Make POT file and no printout on device #7 

1H= 1 

Calculate roll viscous damping 

11= 1 

Non-dimensional output 

1J = 0 

Do not do calculations in wave direction 

1K= 1 

Generate RAO file 

1L = 0 

No shift input offset axis 

1N = 0 

No empirical surge damping values do be read in. 

NS= 13 

Number of input stations 

BPL = 520 

Ship length (ft) 

RO= 1.9905 

Water density (slugs / ft) 

GRAY = 32.174 

Acceleration of gravity (ft / s ) 

CDISPL = 0.0 

Auto check of displacement computed value 

DEPTH = 0.0 

Infinite water depth 

XBKF = 33.0 

x-location of forward end of bilge keel (ft) (not used) 

XBKA = -26.0 

x-location of aft end of bilge keel (ft) (not used) 

BKWTH= 1.0 

Bilge keel width (ft) (not used) 


Table 7. Parameter Settings of SHIPMO.INT'iXq (Part I). 
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Resuming the input settings, the next data is concerned with the ship’s 
body plan. Due to lack of real data, the body plan of a similar ship was used in this simu¬ 
lation, obtained from [10] and depicted in Appendix B. The main structural dimensions 
as length, breadth, draft and average displacement are very much the same as the actual 
ship used in the maneuvering model. The skull was segmented in 13 stations, numbered 
as station 0, 14 , 1, 2, 3, 4, 5, 6, 7, 8, 9, 9 14 and 10. For each station, seven points define a 
curve in the Y-Z plane. This data is used in the two-dimensional potential calculations. 
The two-dimensional potential, y/k, is determined by solving the boundary value problem 
listed in the equation V^^^(y,z) = 0 in the fluid domain, where k may be 1, 2, 3 or 4 cor¬ 
responding to problems for surge, sway, heave or roll, respectively. In this program, the 
boundary value problem is solved by a source distribution technique [9]. The points used 
to represent the hull stations are shown in Fig. 16. The points start at keel and work up. 


Container Ship Body Plan 



•40 -24 -16 -8 -4 0 4 8 16 24 40 

Y(ft) 

Figure 16. Skull Station Points. 


These points are grouped by station sequentially in the input file. For each 
group, the number of the station, its x-position, the number of points used to describe the 
curve and the ILID parameter are posted. The parameter ILID = 0 is used to prevent ir¬ 
regular frequencies. The next parameters to be set are ZCG (vertical location of ship cen- 
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ter of gravity positive above waterline) and RADGRO (roll radius of gyration about x- 
axis). They are set to 9999 and 0.0, respeetively, to indicate that the weight curve com¬ 
puted value will be used. 

The last part of the input file used for this simulation addresses the regular 
wave control. These settings, shown in Table 8, define the number of tables to be gener¬ 
ated in the output fide. In this case, seven speed steps and 13 wave heading angle steps are 
used, generating 7x13 = 91 tables. These choices are directly related to the desired level 
of accuracy in the interpolation process. The smaller the increments, the higher the accu¬ 
racy, and more tables are created. As long as memory space is not a current concern, de¬ 
pending on the application, more or fewer tables may be created. In other words, if accu¬ 
racy is requested for a given application, it will not imply that the processing time of the 
wave model will increase. It will be exactly the same as in applications that do not re¬ 
quire much accuracy. 


Settings 

Description 

WA= 1.0 

Continue wave amplitude 

SWL= 10.0 

Eowest wavelength (ft) 

BWL = 1460.0 

Biggest wavelength (ft) 

DELWL = 50.0 

Increment in wavelength (ft) 

VMIN = 0.0 

Minimum ship speed (ft/s) 

VMAX = 60.0 

Maximum ship speed (ft/s) 

DELV= 10.0 

Increment in ship speed (ft/s) 

WANGI = 0.0 

Initial wave heading angle (degrees) 

WANGA= 180.0 

Einal wave heading angle (degrees) 

DWANG= 15.0 

Increment in wave heading angle (degrees) 


Table 8. Parameter Settings of SHIPMO.IN ViIq (Part II). 

The SHIPMO.IN file used to generate the tables for this simulation is pre¬ 
sent in Appendix C. 
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The output file ereated by this eomputer program, SHIPMO.OUT, eontains 
information about the ship’s structure, stability and the desired tables, one for each 
(speed/wave heading angle) pair, as a function of wavelength values. Since this is a very 
large file, only the initial information and the two first sets of tables are shown in Appen¬ 
dix D to provide an idea of its format. 

b. Converting to waveModelTables.h File 

An auxiliar program called tableConverter.exe was implemented to con¬ 
vert the tables to a readable format for the C++ simulation program. This program reads 
the SHIPMO.OUT file as input and generates a float double-precision type, two- 
dimensional array variable for each table. Each line of this matrix is associated with a 
wavelength value (from 10 ft to 1460 ft in steps of 50 ft), and the column s are: surge am¬ 
plitude, surge phase, heave amplitude, heave phase, pitch amplitude, pitch phase, sway 
amplitude, sway phase, roll amplitude, roll phase, yaw amplitude and yaw phase, in this 
order. These outputs were named Amplitudek and Phaser where k = surge, heave, pitch, 
sway, roll and yaw. The name of each matrix is automatically created by the program and 
contains the speed and wave-heading angle to which the table refers. For example, the ta¬ 
ble named T_Sp40_Whl35 refers to the ship motion parameters related to speed = 40 ft/s 
and wave heading = 135°. 

At the end of the waveModelTables.h file, a matrix named TablePointer 
with the pointers to the first element of each table just converted is created. All these ma¬ 
trices are used to interpolate values in real time during the execution of the simulation. 

4. Interaction with Other Modules 

The four input parameters are needed to produce the 12 outputs that determine the 
instantaneous ship’s position and orientation due the ocean dynamics. The four input pa¬ 
rameters are ship’s speed, wave heading angle, wavelength of the fundamental wave 
component and the elevation of the ocean at the ship’s CG. The input parameters have 
two distinct sources, the ship’s speed comes from the Maneuvering Model module and 
the others come from the Virtual Environment Model module. 
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The instantaneous ship’s speed is eomputed from the eurrent u and v eomponents 
by using the formula speed = +v^ . 

The three parameters from the Virtual Environment module are obtained using 
funetions provided by Vega’s Applieation Program Interfaee (API). The eurrent wave¬ 
heading is taken by the vgGetProp(ocean_calc, VGOCEAN_WAVE_HEADING) function 
and the current wavelength is taken by applying the wave period value obtained from the 
function vgGetProp(ocean_calc, VGOCEAN WAVE PERIOD) in the formula 
wavelength = Tlgjln , where T„ is the wave period and g the gravity acceleration. 

Here, an integration issue between the wave model and Vega arose, concerning 
the wave elevation at ship’s CG parameter. Vega’s API allows getting current wave ele¬ 
vation values of any point in the virtual sea {X-Y plane) but these values are the summa¬ 
tion of all ten cosines components. It is not able to obtain the elevation value from a spe¬ 
cific component. The wave model needs only the current wave elevation of the funda¬ 
mental wave component. In order to solve this problem, a second ocean was created in 
the virtual scenario, 20 meters below the original one. This second ocean, called 
ocean_calc, has the same characteristics of the original one, called ocean_view, but con¬ 
tains only one cosine component, i.e., the fundamental wave component. Therefore, the 
system is now composed by 2 parallel virtual oceans 20 meters apart, but only the 
ocean_view can be seen in the scenario. The ocean calc is used only to get the parame¬ 
ters for calculations. Figure 17 shows a sample of the relationship between the 
ocean view and the ocean calc wave elevations. 
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Ocean-view and Ocean-calc Wave Elevations 



Figure 17. The ocean_view and the ocean_calc Wave Elevations. 

This solution does not imply an important loss in the system performance since 
the second ocean is not rendered. The major ocean computation time is due to rendering 
and texturing. 

Now that a pure cosine wave in the ocean_calc, has been provided, it is needed to 
determine where in the wave the ship’s CG is, in order to apply the corresponding phase 
angle and obtain the correct wave elevation to compute the right dimensional value for 
each output parameter. It was achieved by taking the difference between two measure¬ 
ments of the wave elevation, in a short distance apart, in the direction of the current 
wave-heading angle. It is similar to computing the derivative of the wave signal to get the 
slope of the point of interest (ship’s CG). Since the maximum wave elevation value is ac¬ 
cessible from the current configuration of the ocean, it is possible to compute how far the 
ship’s CG is from the peak of the cosine wave. In other words, this distance corresponds 
to an offset phase angle (a). Figure 18 shows the instantaneous ship’s CG position in the 
wave. 


46 







The Instantaneous Ship CG's Position 



Figure 18. The Instantaneous Ship CG’s Position. 


The angle a is eomputed using a = eos 


WaveEleV(,Q ^ 
WaveElev^^^^ 


and if the slope is posi¬ 


tive, a = 2n - a. 


Onee a is obtained, the wave elevation related to eaeh output parameter ean be 
calculated as follows: 


WaveEleVj^ = WaveElev^^^ ■ cos (a + Phase,^ ) (5.45) 

where k = surge, heave, pitch, sway, roll and yaw. The final values from the wave model 
for each output parameter are obtained applying the dimensional factors. Equation 5.46 
computes the dimensional outputs for the linear motions: surge, sway and heave. 

output = WaveEleVj^ ■ Amplitude (5.46) 


where k = surge, sway and heave. Equation 5.47 computes the dimensional outputs for 
the rotational motions: 


output = WaveEleVi^ ■ Amplitude/^ ■ WaveNumber (5.47) 

where k = pitch, roll and yaw. 

Now, it is necessary to apply a homogeneous transformation to convert the results 
from the Wave Model Coordinate System to Vega Coordinate System. Eigure 19 shows 
the two frames involved in this operation. 
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Figure 19. Wave Model to Vega Coordinate System Conversion. 

The eonversion is carried out by applying a transformation matrix to every posi¬ 
tion generated by the Wave Model. Now, it is required to rotate -90° by the Zw axis. 
Equation 5.48 shows the homogeneous transformation required. 

Xyl To -1 0l 

yy = I 0 0 ■ y^ = . (5.48) 

Zyj [o 0 ij [z^J [ _ 

The parameters related to rotational motions (pitch, roll and yaw) follow the right hand 
rule in both frames. 

The ship’s motion results generated by the Wave Model module are added to the 
results obtained by the Maneuvering Model module (only the four degrees-of-freedom in 
common). Thus, the wave model results do not influence the internal dynamics of the 
Maneuvering Model. A more coupled integration between the two models was attempted. 
Due to the non-linear behavior of the Maneuvering Model, the final result became unsta¬ 
ble. 

5, Validation of the Wave Model I 

The validation of the results of this model was based on two approaches, analysis 
of the curves (where it is possible to check any important discrepancies between the in¬ 
terpolated outputs of the model and the expected real situation) and the visual check us¬ 
ing the graphic visualization provided by the simulator program. In the first approach, 
four sets of curves were generated that represent the Wave Model outputs in four ship’s 
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speeds: 0, 10.7, 21.6 and 32.5 knots. These speeds eorrespond to the steady state values 
of the four engine order regimes: stop, 14 ahead, ys ahead and full ahead. For all those 
eases, the following parameters were used: sea state = 3, dominant wave length = 360 ft, 
and wave heading varying from 0° to 180° in steps of 45°. Eaeh plot presents the wave 
elevation (dashed line) at the ship’s CG and a speeifie output parameter. The objeetive of 
these plots is to show the amplitudes and phase-shifts between the oeean wave and eaeh 
output under different wave heading angles. 

Figure 20 shows the vertieal plane motion parameters (surge, heave and piteh) at 
speed of 0 knots. 


Surge Motion in speed = 0 



Heave Motion in speed = 0 




Figure 20. Vertieal Plane Motion Output (Speed = 0 Knots). 


Observing the plots in Figure 20, it ean be seen that the surge motion is not sensi¬ 
tive to the following and head seas, and praetieally null in the beam sea eondition. As ex¬ 
pected, the heave motion has a very accentuated behavior in the beam sea and a moderate 
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behavior in the other sea conditions. In the pitch motion, it is seen as an accentuated os¬ 
cillation in the quartering and bow seas, and a very small response in the beam sea condi¬ 
tion, as expected. 

Figure 21 shows similar plots to the horizontal plane motion parameters sway, roll 
and yaw. Notice that this model has response zero in the following and head seas for all 
motions in the horizontal plane. The sway motion is more accentuated in the beam sea 
condition, as expected. The yaw motion has a very small response for all wave-heading 
conditions. All these analytical responses are very reasonable considering the type of the 
ship. The visual feedback from the graphics animation was very acceptable also, even in 
high sea state conditions, represented by sea states four to six. 


Sway Motion in speed = 0 





time (s) 


Figure 21. Horizontal Plane Motion Output (Speed = 0 Knots). 


Figures 22 and 23 show the vertical and horizontal plane motion outputs at 
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Surge Motion in speed = 10.7 knots 





Figure 22. Vertical Plane Motion Output (Speed = 10.7 Knots), 
speed of 10.7 knots, respectively. 


Sway Motion in speed = 10.7 knots 
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Figure 23. 


























































From the plots of Figures 22 and 23, it is seen that the outputs have reasonable 
values, whieh eould be attested by the visual eheek. Also, notiee that the wave elevation 
signal inereases in frequeney as the wave heading goes to 180 deg beeause of the Doppler 
effeet. 

Figures 24 and 25 show the vertieal and horizontal plane motion outputs at speed 
of 21.6 knots, respeetively. 


Surge Motion in speed = 21.6 knots 




time (s) 


Figure 24. Vertieal Plane Motion Output (Speed = 21.6 Knots). 


As the speed inereases, an unexpeeted behavior started to oeeur, espeeially in the 
surge and sway motions. It was verified from SHIPMO.OUT file tables that, possibly due 
to some numerieal issues, the amplitude values of some parameters beeame unstable for 
speed values greater than 10 knots. It was also notieed in the visual eheek that, even 
though the plots seem to be eorreet, the ship did not behave as expeeted, espeeially in 
high seas eonditions. 
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Sway Motion in speed = 21.6 knots 
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Figure 25. Horizontal Plane Motion Output (Speed = 21.6 Knots). 

The plots of Figures 26 and 27 show the problem worsening at the maximum 


Surge Motion in speed = 32.5 knots 





Figure 26. Vertical Plane Motion Output (Speed = 32.5 Knots). 
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ship’s speed of 32.5 knots. 


Sway Motion in speed = 32.5 knots 





Figure 27. Horizontal Plane Motion Output (Speed = 32.5 Knots). 

The implemented solution to avoid this kind of instability was to negleet the surge 
and sway eontributions from this wave model when the speed of the ship is greater than 
approximately 20 knots. 

D. THE WAVE MODEL II 

1, Description 

The Wave Model II is a wind generated wave model eoupled to the Maneuvering 
Model. It generates an additional moment to the K component, contributing to the effect 
on the roll motion. This is a second order linear wave model. This kind of approximation 
is usually preferred by ship control systems engineers due to its simplicity and applicabil¬ 
ity. The first applications were reported by Balchen, Jenssen and Sac lid [17]. The 
model’s transfer function is shown below: 

-F (5.49) 

S + “f COq 


54 

























where the gain K„ is defined as = h^co^cr^. The parameter C is the damping eoeffi- 

eient, coo is the dominating wave frequeney and (7„ is the wave intensity eoeffieient. The 
values used for C and ctw were 0.00005 and 1.0, respeetively. 

The model in the state-spaee representation is defined by: 



0 1 

~cOq —2,(^ci)q 



+ 


K 


wit) 


(5.50) 


where w(t) is a zero mean white noise proeess. It introduees a random faetor to the ship’s 
roll motion behavior. The state variable Xhi is the additional moment K'„ave added to Eq. 
5.24 aeeording to the prineiple of superposition. 


2. Interaction with Other Modules 

This model works like a sub-module of the Maneuvering Model module as long 
as they are eoupled. It uses the aetual value of the wave frequeney from the Virtual Envi¬ 
ronment module. The white noise proeess input is internally generated by the mathemati- 
eal funetion randnormalQ. The output parameter K'^ave is used by the Maneuvering 
Model module during the ealeulations of the forees and moments. The new states Xhi and 
Xh 2 are obtained together with the new states of the Maneuvering Model module as in Eq. 
5.43. 


3. Validation of the Wave Model II 

The validation of the Wave Model II was made by repeating the same test used 
when the Maneuvering Model was validated. The test ran in sea state = 3 and the meas¬ 
urements started with the ship at the speed of 21.6 knots. 
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Roll Angle without Wave Model II 




Figure 28. Roll Angle with and without Wave Model II. 

The following maneuvers were made: a 20° rudder angle turn to port followed by 
a -20° rudder angle turn to starboard, and then returning to the neutral rudder position. 
Figure 28 shows the two test eonditions. It is elear the effeet of the prineiple of superposi¬ 
tion in the roll motion when the Wave Model II is eoupled to the Maneuvering Model. 
The use of the Wave Model II ereated a random behavior in the ship motion, making the 
final result more realistie. The introduetion of the K'„ave component eaused a small dis- 
turbanee in the Maneuvering Model’s yaw motion. However, this disturbanee does not 
produee any unstable behavior in the model and it is not pereeptible in the visual eheek. 

E. THE WIND MODEL 

1. Description 

The resultant wind forees and moments aeting on a surfaee vessel are usually de¬ 
fined in terms of relative wind speed Vr (knots) and relative angle jr (deg) that are de¬ 
fined by: 



(5.51) 

( \ 


-1 '^R 

7 ^ = tan ^ 

(5.52) 
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where the eomponents of Vr in the x andy direetions are defined by: 


+ (5.53) 

VR=V^^m{y)-v + v^. (5.54) 

Here, u and v are the ship veloeity eomponents, the Uc and are the oeean eurrent 
veloeity eomponents, ¥„ is the wind speed while y = (//»- i// is the angle of relative wind 
of the ship bow. 

Typieally wind models only treat the foree and moments that are direetly related 
to surge, sway and yaw motions, i.e., = . The model adopted in 

this study is one suggested by Isherwood [18], whieh eonsiders wind resistanee of mer- 
ehant ships. This model defines the wind forees for surge and sway, in Newtons (N), and 
the wind moment for yaw, in Newton-meters (Nm) as is shown in Eqs. 5.55-57, respee- 
tively: 

(5.55) 

^.^XcAr,)p.rlA, (5.56) 

N,,^=\c,{r,)pJlA,L (5.57) 

where Cx, Cj and Cn are the foree and moment eoeffieients, pw is the density of the air in 
kg/m , At and Al are the transverse and lateral projeeted areas in m and L is the overall 
length of the ship in meters. Based on these equations, measured data were analyzed by 
multiple regression teehniques in terms of eight parameters: length (T), beam {B), lateral 
projeeted area {Al), transverse projeeted area {At), lateral projeeted area of superstrueture 
{Ass), length of perimeter of lateral projeetion of model exeluding waterline and slender 
bodies sueh as masts and ventilators (5), distanee from bow of eentroid of lateral pro¬ 
jeeted area (Q, and number of distinet groups of masts or kingposts seen in lateral pro¬ 
jeetion {Af). Furthermore, Isherwood [18] found that the data were best fitted to the fol¬ 
lowing equations: 


57 





- A + 4 


24 24 L S . C . 

"^ 2 -^ +4 —+ ^4 —+ 4 Y"*" "^6^ 


(5.58) 


^ „ r.2A, ^ 2A^ ^ L ^ S „ C ^ A,, 

Cv — B^ + B, ^ 1- B-, ;;—l- B^ —l- B, —l- B^ —h B^ 


L 


2 ' ^2 ^2 —3 


B 


L 


L 


A, 


(5.59) 


c„ = c+c^+a^+c-+c-+c^ 
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(5.60) 


where the coefficients Ai and Bi (/ = 0.. .6) and Q (/ = 0.. .5) come from tables in function 
of the relative angle jr, that varies from 0 to 180 degrees. These tables were implemented 
in the simulation program and all coefficients are computed in real time by first order in¬ 
terpolation according to the current relative angle jr. 

The values used for the eight parameters related to the container ship of this study 
are shown in the Table 9. 


Parameter 

Value 

L 

175 m 

B 

25.40 m 

Al 

1400.0 

At 

203.3 m^ 

Ass 

200.0 m^ 

S 

16 m 

c 

90 m 

M 

1 


Table 9. Ship Parameters for the Wind Model. 


In order to be used correctly, the parameters Y„ind and Nwind were converted 
to their non-dimensional forms, Ywind and N'mnd, by applying the associated factors 
described in Table 3, i.e., the force parameters were divided by pU L /2 and the moment 

2 3 

parameter was divided hy pU L /2. 

2 . Interaction with Other Modules 

The wind model works like a sub-module of the Maneuvering Model module as 
long as they are coupled. The input parameters come from the Environment Model mod- 
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ule. They are the u and v eomponents of the relative wind, and the eurrent relative angle 
Jr. The Maneuvering Model module uses the output parameters X'm„d, Ywind and 
during the ealeulations of the forees and moments. 


3, Validation of the Wind Model 

The validation of the Wind Model was made through the observation of the ship 
behavior under several situations. One test is reported here to allow its analysis. The ini¬ 
tial eonditions were the ship’s heading angle = 0° and the wind’s heading angle = 270°, 
i.e., the wind is eoming from starboard. The wind intensity was set as 50 knots in order to 
make the measurements easier in terms of magnitude of the values and the duration of the 
test. It was observed that the ship tends to be aligned with the wind, i.e., the wind makes 
the ship’s bow move to the opposite direetion in the wind. This movement is seen from 
Figure 29 in the roll and yaw plot. The relative angle yr starts at nil rad (wind eoming 
from starboard) and goes to % rad (head wind). The ship’s heading angle starts at 0 rad 
and goes to nil rad (90°), pointing exaetly to where the wind is eoming from, at the end 
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Figure 29. Wind Model Test. 


of the simulation. Sinee the models are eoupled, a roll movement was indueed too, as 
seen in the same plot. The amplitude of the roll motion inereased as the simulation time 

59 













increased, but when the ship beeame aligned with the wind, the roll amplitude started to 
deerease, as expeeted. From the veloeities plot in the same figure, notiee that the ship is 
rotating until about 18 minutes of simulation and then the ship is pushed baekwards by 
the wind. An oscillated behavior is seen in the v veloeity component. It is associated to 
the induced roll motion. 

The observed behavior of the ship during the simulation was as expected and the 
model was validated. 


F. THE OCEAN CURRENT MODEL 
1. Description 

The ocean eurrent model used in this study is a simple two-dimensional eurrent 
model. It is eommonly used in surfaee vessel applications. This model considers the loeal 
effect of the current, applied to the ship’s CG. Therefore, this model has two parameters 
only, which are the average eurrent speed Vc and direction of the current The two ve¬ 
locity components in the Earth-fixed frame can be written as: 

=V^cos/3 (5.61) 

v^=V^smj3. (5.62) 


Considering the horizontal motion only, those eomponents ean easily be con¬ 
verted to the fixed-body frame using the transformation matrix below: 


V„ 


eos^ 

sin^ 


-sin^ 

cost// 

_vf_ 


(5.63) 


where y/ is the heading angle of the ship. Substituting Eqs. 5.61 and 5.62 into Eq. 5.63 fi¬ 
nally yields: 

cos(/3-y/) (5.64) 

sm(/3-y/). (5.65) 


The parameters Vc and ^ ean be ehanged by the user at any time during the simu¬ 
lation. 
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2, Interaction with Other Modules 

The ocean current heading angle and velocity parameters are obtained from the 
Environment Model module. The outputs of the model, the velocity components Uc and 
Vc, are then integrated as in Eq. 5.43, to produce the new x andy position values after the 
integration interval of time At. The velocity components Uc and Vc also are used to set the 
input parameters to the Wind Model module as shown in Eqs. 5.53 and 5.54. 

The ocean current has a graphic representation in the simulation scenario. The 
two buoys are presented at the beginning of the simulation, located on both sides of the 
ship. Their eddies are rendered as a function of the parameters Vc and p. Ear from the ini¬ 
tial point of the simulation, the effects of the ocean current are more easily noticed if the 
ship is just drifting. Otherwise, the user can check the parameter values at any time on the 
screen in text format. 

3, Validation of the Ocean Current Model 

The validation of this model was very simple. Erom the initial point of simulation, 
with the ship in the drifting condition, the ship motion can be easily verified in the direc¬ 
tion of the ocean current heading angle, related to the two fixed buoys. 

This chapter presented a theoretical discussion about hydrodynamic models and 
described in detail all the motion models used in this thesis. The following chapter is re¬ 
lated to the software implementation. The simulation system is divided into blocks and 
the software developed in each block is described. 
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VI. THE SIMULATION PROGRAM 


This chapter shows the structure of the simulation program. All program eompo- 
nents and the relationship between them are described. 

A, GENERAL DESCRIPTION 

The main purpose of this simulation program was to provide the user a virtual 
laboratory to test the maneuvering of a marine vessel under a pre-defined set of environ¬ 
ment disturbances through the observation of its behavior. The program provides re¬ 
sources to help the user during the simulation. Three arrows that represent eaeh direction 
vector for wind, ocean current and wave heading, a 0.1 nautical-mile grid, a compass and 
a stopwatch can be activated by the user at any time during the simulation. 

The program was built in the Mierosoft Visual Studio 6.0 environment. It is a 
multi-thread C++ application and the threads communicate via shared memory. The 
simulation program runs as an application layer over the Vega system, using Vega’s API 
(application program interface). Figure 30 shows the structure of the simulation program. 



Screen Commands 

Figure 30. Simulation Program Structure. 
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The Main Program block creates the threads and the events needed to control the 
access to the shared memory area. In this thesis, only one vehicle (container ship) is im¬ 
plemented, but the program could accept the introduction of more vehicles with few 
changes. 

The Simulation Controller block {simulator_controller thread) is responsible for 
controlling all the operation of the simulator. Only this thread interfaces with Vega’s 
API. Therefore, this block controls the graphic presentation, the input of user commands 
and the environment management. 

The Container Ship block {container_shipjnotion thread) contains all five model 
implementations. The following sections describe in details these two blocks. 


B, THE SIMULATION CONTROLLER BLOCK 

The content of the Simulation Controller block is depicted in Figure 31. Now, the 
principal functions of each one of its elements are described below. 
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Figure 31. Simulation Controller Scheme. 
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1. Initialization 

This block does the initialization of the simulation program. The initialization is 
divided in tasks as follows. At first, Vega and its speeial modules (Vega Marine and 
Vega Speeial Effeets) have to be initiated. After that, the containers hip _simulator.ADF 
file has to be loaded to inform Vega about all the pre-made work exeeuted in Lynx. As it 
was mentioned before. Lynx is a graphieal interfaee that allows the user to eonstruet the 
bulk of the applieation in a point and eliek environment. Lynx generates an Applieation 
Definition Lile (ADL) that is used to eonfigure the applieation. Onee the file is loaded, 
the simulation program is able to obtain all the initial settings by using the Vega’s API 
funetion ealls, to initiate the internal variables in order to use during the program exeeu- 
tion. 


2. Main Loop 

This bloek is an infinite loop that lasts until the simulation exists. It is eomposed 
of five distinet phases: 

• Data aequisition from the virtual environment in order to prepare the 
model’s input parameters (Data Aequisition bloek); 

• Information exehange through the eritieal seetion or shared memory. Here, 
the model input parameters are written and the model output results from 
the last iteration are read (Shared Memory Transfer I bloek); 

• User’s eommands and settings (Input New Settings bloek); 

• Position and orientation updating for the objeets of the seene (Positions 
and Orientations bloek); 

• Rendering permission eommand to Vega (Rendering Control bloek). 

3. Data Acquisition 

This bloek reads variables of interest from the eurrent state of the environment 
and prepares the model inputs. There are three phases in this operation. The variables ob¬ 
tained in eaeh phase are shown below: 

• Wave parameters: wave_period (period of the fundamental wave eompo- 
nent), wave_heading_R (wave heading related to the ship heading), and 
alpha (phase angle of the ship CG position related to the wave_elev_max 
point). All oeean readings are referred to the ocean_calc; 
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• Ocean current parameters: curr_heading_R (current heading related to the 
ship heading), u current R and v_current_R (u and v components of the 
current velocity in the Maneuvering Model coordinate frame); 

• Wind parameters: wind_gamma_R (wind heading related to the ship head¬ 
ing), u_current_R and v_current_R {u and v components of the wind ve¬ 
locity in the Maneuvering Model coordinate frame). 

4, Shared Memory Transfer I 

The following variables are written on the shared memory: ship_status 
(DRIFTING, CRUISING or SINKING), framejime (integration interval), shaft_order 
and rudder_order, wave_period, wave_heading_R, alpha, wave_elev_max, u current R 
and v_current_R, wind_gamma_R, u current_R and v_current_R, aux_l, aux_2 and 
aux_3 (auxiliary variables for ploting), and trigger (starts/finishes the recording of the log 
fde). 

The following variables are read from the shared memory: ship_x, ship_y and 
ship_z (positions), ship_h, ship_r, ship_p (orientations), ship_u and ship_v (velocities), 
ship_n (actual shaft speed), ship_delta (actual rudder angle), and speed (actual ship 
speed). 

After that, the container_shipjnotion thread is waked to provide the next set of 
model outputs. 


5, Input New Settings 

This block is responsible for the input of commands and the consequent new set¬ 
tings in the Vega context, as required. The commands are separated in groups as follows: 

• Visualization commands: change the current point of view into perspec¬ 
tive, port side, astern, bow, bridge and above positions, and execute zoom- 
in and zoom-out in each of them (except bridge); 

• Ship commands: change the engine and rudder orders; 

• Information commands: toggle the grid, the compass, the vectors and the 
text on the screen, and starts/stops/clears the stopwatch; 

• Environment commands: set the sea state, toggle the fog, increase and de¬ 
crease the daylight, the intensity of the wind and the ocean current, and 
change the direction of the wind, the ocean current, and the wave-heading; 
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• Auxiliary commands: increase/decrease values of test variables, start/stop 

the reeording of the log file and display variables of interest in the system 

eommand window. 

6. Positions and Orientations 

This bloek implements the ship status eheeking. The ship status eould be 
DRIFTING, CRUISING or SINKING. Based on the ship status information, this bloek 
updates the position and orientation of the ship and the direetion of the smoke. The posi¬ 
tion and orientation of other objeets of the seene as the eompass and the direetion veetors 
are also updated here. 

7. Rendering Control 

The Rendering Control bloek synehronizes frame proeessing using the Vega’s 
API vgSyncFrame(), and starts the rendering of a new frame with vgFrame(). In this 
bloek, all variables that will show up on the sereen must be updated in order to be used 
by the post draw drawHUD() funetion. Here, if the stopwateh is in the COUNTING state, 
it is also updated based on the frame time. 

C. THE CONTAINER SHIP BLOCK 

The Container Ship bloek {container_ship_motion thread) contains all five model 
implementations as shown in Figure 32. 

1. Initialization 

This bloek does the initialization of the eontainer ship motion models. All state 
variables are initialized here. 

2. Main Loop 

This bloek is an infinite loop that lasts until the simulation exists. It is composed 
of four distinet phases: 

• Information exehange through the eritieal seetion or shared memory. Here, 
the model input parameters are read and the model output results from the 
last iteration are written (Shared Memory Transfer II bloek); 

• Exeeution of the model eode; 
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Numerical integration; 
Coordinate frame conversion. 







Figure 32. Container Ship Block. 

3. Shared Memory Transfer II 

The following variables are read from the shared memory: ship_status 
(DRIFTING, CRUISING or SINKING), framejime (integration interval), shaft_order 
and rudder_order, wave_period, wave_heading_R, alpha, wave_elev_max, u current R 
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and v_current_R, wind_gamma_R, u current R and v_current_R, aux_l, aux_2 and 
aux_3 (auxiliary variables for ploting), and trigger (starts/fmishes the reeording of the log 
fde). 

The following variables are written in the shared memory: ship_x, shipjy and 
ship_z (positions), ship_h, ship_r, ship_p (orientations), ship u and ship_v (veloeities), 
ship_n (aetual shaft speed), ship_delta (aetual rudder angle), and speed (aetual ship 
speed). 

After that, the simulation_controiler thread is waked up to update the graphieal 
seene with the last eomputed data. 

4, Coupled Models 

This bloek holds three eoupled models: the Maneuvering Model, the Wind Model 
and the Wave Model II. It is represented by the containerShipCoupledModels (x, x_dot, 
input) funetion. This funetion has the state veetor x and the veetor input as inputs and the 
state derivative veetor x_dot as output. It is ealled in the drifting and eruising ship state 
situations. The Wind Model uses the set of tables deelared in the windModelTables.h file. 
Those tables are interpolated to eompute the Aq.. 6 , Bo ,,6 and Co ..5 coeffieients. 

5, Ocean Current Model 

This bloek simply gets the surge and sway veloeity eomponents and passes them 
to the Numerieal Integration bloek. 

6, Wave Model I 

This bloek is represented by the waveModel (inputWM, outputWM) funetion. It is 
ealled in all ship state situations (drifting, eruising and sinking). The Wave Model I bloek 
uses the set of tables deelared in the waveModelTables.h file. 

7, Numerical Integrator 

This block implements the Euler’s integration method to compute the next state 
variable values for the Coupled Models block and the Ocean Current Model block. The 
Numerical Integrator block uses the framejime variable as integration interval. 
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8, Coordinate Frame Converter 

This block implements the required eoordinate transformations to adjust the 
model outputs to the eoordinate frame used by the Vega system. Here, the model outputs 
are eombined in order to produee the 6-DOF paek of variables that represent the position 
and orientation of the ship in the next graphieal frame. 

This ehapter deseribed how the simulation system software is organized. The next 
ehapter shows how to operate the simulation program, deseribing all the eommands the 
operator ean use. 
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VII. OPERATION OF THE SIMULATOR 


This chapter describes the operational procedures of the ship simulation system 
developed in the previous chapters. 

A, GETTING STARTED 

The simulation program is started by running the shipSimulator.exe executable 
file. The simulation window opens, showing the container ship drifting between two 
buoys. The engine is running (indicated by the ship smoke), but the clutch is not into 
gear. The following initial conditions are observed: 

• Point of view: perspective; 

• Fog: ON and daylight: maximum; 

• Sea state: 1 and wave heading: 135 degrees; 

• Ocean current intensity: 0 knots and direction: 0 degrees; 

• Wind intensity: 0 knots and direction: 225 degrees; 

• Grid: OFF, direction vectors: OFF, and compass: OFF; 

• Stopwatch: reset; 

• Screen text: ON. 

B, COMMANDS 

The operator introduces commands using the keyboard. They are grouped in four 
sets: visualization, ship order, information, environment and auxiliary commands. Table 
10 shows the visualization commands. 


Key 

Action 

P 

Changes the point of view position to perspective, port 
side, astern, bow, bridge and above positions. 

z 

Increase zoom out 

Z 

Increase zoom in 


Table 10. Visualization Commands. 
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Table 11 shows the ship order commands. 


Key 

Action 

Up arrow 

1 2 

Increase engine order from Stop, f Ahead, f 
Ahead to Full Ahead 

Down arrow 

Decrease engine order from Full Ahead, f Ahead, 
V 3 Ahead to Stop 

Right arrow 

Turn rudder to starboard up to 20 deg (1 deg steps) 

Left arrow 

Turn rudder to port up to -20 deg (1 deg steps) 


Table 11. Ship Order Commands. 


Table 12 shows the environment commands. 


Key 

Action 

/ 

Toggle the fog 

/ 

Decrease the daylight 

L 

Increase the daylight 

W 

Increase the wind direction (5 deg steps) 

w 

Decrease the wind direction (5 deg steps) 

C 

Increase the ocean current direction (5 deg steps) 

c 

Decrease the ocean current direction (5 deg steps) 

H 

Increase the wave heading angle (5 deg steps) 

h 

Decrease the wave heading angle (5 deg steps) 

E 

Increase the wind intensity up to 50 knots (5 knots steps) 

e 

Decrease the wind intensity down to 0 (5 knots steps) 

V 

Increase the current intensity up to 3 knots (0.5 knots steps) 

V 

Decrease the current intensity down to 0 (0.5 knots steps) 

1 

Set the sea state to 1 

2 

Set the sea state to 2 

2 

Set the sea state to 3 

4 

Set the sea state to 4 

5 

Set the sea state to 5 

6 

Set the sea state to 6 


Table 12. Environment Commands. 
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Table 13 shows the information commands. 


Key 

Action 

g 

Toggle the grid 

a 

Toggle the direction vectors 

q 

Toggle the compass 

t 

Toggle the text on the screen 

s 

Command the stopwatch 


Table 13. Information Commands. 


The auxiliary commands are used for testing, plotting and debugging purposes. 
They are used during the development phase to make easier to conduct visual tests, to 
create a log file for plotting and to help in debugging operations. Table 14 shows the aux¬ 
iliary commands. 


Key 

Action 

M/N 

Increase the value of a variable to check its visual effect 

m/n 

Decrease the value of a variable to check its visual effect 

r 

Start/stop the recording of the log file for plotting purpose 

d 

Display variables of interest on the system’s window 


Table 14. Auxiliary Commands. 


This chapter described all the commands required to operate the simulation pro¬ 
gram properly. The next chapter ends this study with the thesis conclusions and sugges¬ 
tions for future work. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 


The objective of this thesis was to develop a computer simulation program in 
which physical models were implemented in order to achieve a realistic representation of 
a ship in a virtual environment considering its physical features such as the engine, pro¬ 
peller and skull in the presence of environment conditions including waves, ocean cur¬ 
rent, wind, fog and day/night issues. This objective was reached. A simulation program 
was developed and it worked well. All the modeling goals defined by the initial survey 
were accomplished, except the weight distribution response. However, this item was of 
less relevance in the expert’s point of view for naval tactical simulators. The program is a 
good starting point for further implementations 

Complete marine models are difficult to find in the literature probably because 
their construction depends on huge experimental efforts to define all the hydrodynamic 
derivative coefficients. However, this study showed that the integration of multiple mod¬ 
els from different sources is feasible and produces good results, meeting the application 
requirements. 

The use of the interpolation technique to take advantage of complex models 
worked well. The linear interpolation of pre-assembled numerical tables produces a very 
low computational load when running in real-time and the results were very acceptable. 

The Vega system showed to be very useful to build complex virtual environment 
applications in a short time. Some special effects from Vega Marine module are very 
good and the ocean simulation is fine. The bow wave effect, however, does not work well 
especially in high seas conditions. The number of API functions is large and they are easy 
to use. The built-in manual is very useful too. 

The use of dynamic models with relative low complexity produced sufficiently 
good results, saving processing time for other tasks in complex simulation systems. 

As a recommendation for future work, a natural way of continuing this thesis 
work would be the integration with sensor models, such as radar models, that consider the 
sea state and the weather conditions as inputs to impose some constrains in their usage. 
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Weapons could be included in the experiment as well. Naval tactical simulators that in¬ 
clude all these features should be very useful to make the officer’s training more realistic 
and immersive. 
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APPENDIX A. SURVEY 


If you were going to use a Simulator for Naval Tactieal Training in a 3D Virtual 
Environment with the following features: 

• 3D graphical scenery of operations where you are able to visually track the 
maneuvers of the participant platforms (ships, helicopters, airplanes, sub¬ 
marines) from several points of view. 

• All platforms are equipped with full armament and detection equipments. 

• The system allows engagement between platforms. 

Question ; What would you consider IMPORTANT MODELING 
PARAMETERS in terms of physical characteristics of the SHIPS 
and the ENVIRONMENT that could make the TRAINING more re¬ 
alistic? 

Note 1 ; Please, rank the options below starting from the most to the least relevant 
in your opinion, 1 being most relevant. More than one option can have 
the same level of relevance. 

Note 2 ; In your evaluation, consider PRIORITIZING the enhancement of 
FUNCTIONAL issues of the Simulator over esthetic issues. 

_The ship behavior in terms of acceleration/deceleration rates 

_The ship behavior in terms of turning rate 

_The ship behavior in terms of weight distribution (tilt/trim) 

_The wake of the ship 

_The influence of sea state (ship “rocking”) 

_The influence of currents in the ship behavior 

_The influence of the wind in the ship behavior 

_Visibility issues (day/night, fog) 

_Other:_ 

Your area of activity ;_ 
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PRINCIPAL DIMENSIONS 


APPENDIX B. CONTAINER SHIP BODY PLAN [AFTER REF. 10] 
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HALF-BREADTH PLAN 
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APPENDIX C. SHIPMO.IN INPUT FILE 


The SHIPMO.IN file showed below is the input file to generate the numerieal ta¬ 
bles by running the program. 


Container ship model 

0000101110100 13 0 

520.0000 1.9905 32.1740 1.26E-05 0.0 0.0000 

33.0000 -26.0000 1.0000 

7 260.0000 0.0000 0 

0.0000 -29.0000 

1.1200 -28.5000 

2.5300 -25.0000 

2.4900 -21.0000 

1.1200 -13.0000 

0.0100 - 2.0000 
0.5600 0.0000 

7 234.0000 0.0000 0 

0.0000 -29.0000 

1.9600 -28.5000 

4.2100 -25.0000 

4.7700 -21.0000 

4.2100 -13.0000 

4.4900 -2.0000 

4.8000 0.0000 

7 208.0000 0.0000 0 

0.0000 -29.0000 

3.3700 -28.5000 

6.7400 -25.0000 

8.1400 -21.0000 

8.4200 -13.0000 

9.8200 -2.0000 

10.1600 0.0000 

7 156.0000 0.0000 0 

0.0000 -29.0000 

6.1800 -28.5000 

13.1900 -25.0000 

16.8400 -21.0000 

19.6500 -13.0000 

23.0200 -2.0000 

24.0000 0.0000 

7 104.0000 0.0000 0 

0.0000 -29.0000 

12.9100 -28.5000 

24.0000 -25.0000 

27.2300 -21.0000 

30.3200 -13.0000 

33.1200 -2.0000 

33.3200 0.0000 

7 52.0000 0.0000 0 

0.0000 -29.0000 

24.9800 -28.5000 

33.1200 -25.0000 

35.3700 -21.0000 

36.7700 -13.0000 

37.8900 -2.0000 

37.8900 0.0000 

7 0.0000 0.0000 0 

0.0000 -29.0000 

29.4700 -28.5000 

35.6500 -25.0000 

37.0500 -21.0000 
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.8900 

-13 

37 . 

.8900 

-2 
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.8900 

0 

7 

- 52 . 

.0000 

0 . 

.0000 

-29 

24 . 

.9800 

-28 

33 . 

.1200 

-25 

35 . 

.3700 

-21 

36 . 

.7700 

-13 

37 . 

.8900 

-2 

37 . 

.8900 

0 

7 

-104 . 

.0000 

0 . 

.0000 

-29 

18 . 

.9200 

-28 

26 . 

.5400 

-25 

30 . 

.4900 

-21 

35 . 

.0100 

-13 

37 . 

.2700 

-2 

37 . 

.5500 

0 

7 

- 156 . 

.0000 

0 . 

.0000 

-29 

6 . 

.7800 

-28 

14 . 

.4000 

-25 

19 . 

.2000 

-21 

25 . 

.4100 

-13 

32 . 

.4700 

-2 

33 . 

.6000 

0 

7 

- 208 . 

.0000 

0 . 

.0000 

-29 

2 . 

.2500 

-28 

5 . 

.3600 

-25 

7 . 

.3400 

-21 

11 . 

.0100 

-13 

20 . 

.0500 

-2 

22 . 

.0200 

0 

7 

-234 . 

.0000 

0 . 

.0000 

-29 

1 . 

.1300 

-28 

2 . 

.2600 

-25 

2 . 

.8200 

-21 

3 . 

.6700 

-13 

11 . 

.8600 

-2 

13 . 

.8400 
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2 - 260 , 

,0000 

0.0000 

-2 

4.8000 

0 
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0 
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10 
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0000 

0000 
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0000 
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0000 

0000 

0000 

0000 

0000 

0.0000 0 

0000 
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0000 
0000 
0000 
0000 
0000 

0.0000 0 

0000 
5000 
0000 
0000 
0000 
0000 
0000 

0.0000 0 

0000 
5000 
0000 
0000 
0000 
0000 
0000 

0.0000 0 

0000 
5000 
0000 
0000 
0000 
0000 
0000 


0.0000 0 


0000 


0000 


0000 


o 

o 

o 

o 

0000 

1500.0000 

0000 

15.0000 


50.0000 00 


0000 60.0000 


10.0000 
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APPENDIX D. SHIPMO.OUT OUTPUT FILE 


The SHIPMO.OUT file is the output file of the SHIPMO.BM Sinee it is 

a very large file, only the initial information and the two first sets of tables (of a total of 
91 tables) are shown below to provide an idea of its format. 


IContainer ship model 

OOPTION CONTROL TAGS - ABCDEFGHIJKLMNNS NPAC 

000010111010 13 0 13 0 

0***BASIC INPUT DATA*** 

******* INFINITE DEPTH ******* 

0 LENGTH= 520.000 

DENSITY= 1.9905 GAMMA= 64.0423 GNU= 0.126000E-04 GRAVITY= 32.1740 
OBILGE KEEL DATA; XFOR= 33.0000 XAFT= -26.0000 LENGTH= 59.0000 WIDTH= 1.0000 


OSTA NO 

XAXIS 

1/2BEAM 

DRAFT 

AREA 

AREA COEF 

ZBAR 

BILGE R 

AREA COEF2 

1 

260.0000 

0.5600 

29.0000 

75.8650 

2.3357 

-18.9654 

0.0000 

2.3357 

2 

234.0000 

4.8000 

29.0000 

244.6150 

0.8786 

-13.8322 

0.0000 

0.8786 

3 

208.0000 

10.1600 

29.0000 

469.6700 

0.7970 

-13.0276 

0.0000 

0.7970 

4 

156.0000 

24.0000 

29.0000 

1046.3350 

0.7517 

-12.5382 

0.0000 

0.7517 

5 

104.0000 

33.3200 

29.0000 

1631.6801 

0.8443 

-13.1316 

0.0000 

0.8443 

6 

52.0000 

37.8900 

29.0000 

2039.7400 

0.9282 

-13.8060 

0.0000 

0.9282 

7 

0.0000 

37.8900 

29.0000 

2118.1150 

0.9638 

-14.0585 

0.0000 

0.9638 

8 

-52.0000 

37.8900 

29.0000 

2039.7400 

0.9282 

-13.8060 

0.0000 

0.9282 

9 

-104.0000 

37.5500 

29.0000 

1865.4099 

0.8565 

-13.2276 

0.0000 

0.8565 

10 

-156.0000 

33.6000 

29.0000 

1337.6200 

0.6864 

-11.8926 

0.0000 

0.6864 

11 

-208.0000 

22.0200 

29.0000 

651.1600 

0.5098 

-10.4750 

0.0000 

0.5098 

12 

-234.0000 

13.8400 

29.0000 

306.9000 

0.3823 

-9.2436 

0.0000 

0.3823 

13 

-260.0000 

4.8000 

2.0000 

9.6000 

0.5000 

-0.6667 

0.0000 

0.5000 


OTHE FOLLOWING STATIONS AND OFFSETS ARE USED FOR TWO-DIMENSIONAL CALCULATIONS 
STA NO= 1 XAXIS= 260.0000 NT= 7 ILID= 0 
0.0000 -29.0000 

1.1200 -28.5000 

2.5300 -25.0000 

2.4900 -21.0000 

1.1200 -13.0000 

0.0100 - 2.0000 
0.5600 0.0000 

STA NO= 2 XAXIS= 234.0000 NT= 7 ILID= 0 
0.0000 -29.0000 

1.9600 -28.5000 

4.2100 -25.0000 

4.7700 -21.0000 

4.2100 -13.0000 

4.4900 -2.0000 

4.8000 0.0000 

STA NO= 3 XAXIS= 208.0000 NT= 7 ILID= 0 
0.0000 -29.0000 

3.3700 -28.5000 

6.7400 -25.0000 

8.1400 -21.0000 

8.4200 -13.0000 

9.8200 -2.0000 

10.1600 0.0000 

STA NO= 4 XAXIS= 156.0000 NT= 7 ILID= 0 
0.0000 -29.0000 

6.1800 -28.5000 

13.1900 -25.0000 

16.8400 -21.0000 

19.6500 -13.0000 

23.0200 -2.0000 

24.0000 0.0000 
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STA N0= 5 XAXIS= 104.0000 NT= 7 ILID= 0 

0.0000 -29.0000 

12.9100 -28.5000 

24.0000 -25.0000 

27.2300 -21.0000 

30.3200 -13.0000 

33.1200 -2.0000 

33.3200 0.0000 

STA NO= 6 XAXIS= 52.0000 NT= 7 ILID= 0 

0.0000 -29.0000 

24.9800 -28.5000 

33.1200 -25.0000 

35.3700 -21.0000 

36.7700 -13.0000 

37.8900 -2.0000 

37.8900 0.0000 

STA NO= 7 XAXIS= 0.0000 NT= 7 ILID= 0 

0.0000 -29.0000 

29.4700 -28.5000 

35.6500 -25.0000 

37.0500 -21.0000 

37.8900 -13.0000 

37.8900 -2.0000 

37.8900 0.0000 

STA NO= 8 XAXIS= -52.0000 NT= 7 ILID= 0 

0.0000 -29.0000 

24.9800 -28.5000 

33.1200 -25.0000 

35.3700 -21.0000 

36.7700 -13.0000 

37.8900 -2.0000 

37.8900 0.0000 

STA NO= 9 XAXIS= -104.0000 NT= 7 ILID= 0 
0.0000 -29.0000 

18.9200 -28.5000 

26.5400 -25.0000 

30.4900 -21.0000 

35.0100 -13.0000 

37.2700 -2.0000 

37.5500 0.0000 

STA NO= 10 XAXIS= -156.0000 NT= 7 ILID= 0 
0.0000 -29.0000 

6.7800 -28.5000 

14.4000 -25.0000 

19.2000 -21.0000 

25.4100 -13.0000 

32.4700 -2.0000 

33.6000 0.0000 

STA NO= 11 XAXIS= -208.0000 NT= 7 ILID= 0 
0.0000 -29.0000 

2.2500 -28.5000 

5.3600 -25.0000 

7.3400 -21.0000 

11.0100 -13.0000 

20.0500 -2.0000 

22.0200 0.0000 

STA NO= 12 XAXIS= -234.0000 NT= 7 ILID= 0 
0.0000 -29.0000 

1.1300 -28.5000 

2.2600 -25.0000 

2.8200 -21.0000 
3.6700 -13.0000 

11.8600 -2.0000 
13.8400 0.0000 

STA NO= 13 XAXIS= -260.0000 NT= 2 ILID= 0 
0.0000 - 2.0000 
4.8000 0.0000 

OCOMPARISON OF INPUT AND COMPUTED DATA 

INPUT DISPL= 19648.7637 XCG= -7.5667 RELATIVE TO ORIGIN AT MIDSHIP 

HYDROSTSTIC DISPL= 19648.7637 XCG= -7.5667 RELATIVE TO ORIGIN AT MIDSHIP 
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OBLOCK COEFFICIENT= 0.6014 
OSTABILITY PARAMETERS 

LCF = -18.1800 ZCG 

WATERPLANE AREA = 

BML = 598.0724 GML 

BMT = 15.5936 GMT 

ORADII OF GYRATION 

KYY = 130.0000 KXX ^ 

IContainer ship model 


-2.0000 ZCB = -13.1960 

28801.2402 
586.8764 
4.3976 

26.5230 PRODUCT OF INERTIA(I46) 


0.207022E+08 


***CONDITIONAL INPUT DATA*** 

WAVE AMPLITUDE= 1.0000 

INITIAL WAVELENGTH= 10.0000 FINAL WAVELENGTH= 1500.0000 DELTA WAVELENGTH= 50.0000 
INITIAL VEL= 0.0000 FINAL VEL= 60.0000 DELTA VEL= 10.0000 
INITIAL WAVE HEADING ANGLE= 0.0000 FINAL WAVE HEADING ANGLE= 180.0000 

DELTA WAVE HEADING ANGLE= 15.0000 


***UNITS OF OUTPUT*** 

NONDIMENSIONAL OUTPUT 
OUTPUT IS DIVIDED BY; 

LINEAR MOTIONS - WAVE AMPLITUDE (WA) 

ROTATIONAL MOTIONS - WAVE SLOPE (WA*WAVEN) 

VELOCITIES - WA*SQRT(GRAV/BPL) 

ACCELERATIONS - WA*GRAV/BPL 
SHEAR - WA*GAMMA*BEAM*BPL 
BENDING MOMENTS - WA*GAMMA*BEAM*BPL**2 
IContainer ship model 

0 SPEED = 0.0000 WAVE ANGLE = 0.00 DEG. 

+ VERTICAL PLANE RESPONSES 

+ (NON-DIMENSIONAL) 


WAVE 

ENCOUNTER WAVE 

WAVE/SHIP 










S U 

R G E 

H E 

A V E 

P I 

T C H 

F R E Q U 

E N C I E 

S LENGTH 

LENGTH 

AMPL, 

PHASE 

AMPL. 

PHASE 

AMPL. 

PHASE 

4.49617 

4.49617 

10.000 

0.0192 

0.0000 

-173.2 

0.0000 

-148.8 

0.0000 

-154.2 

1.83555 

1.83555 

60.000 

0.1154 

0.0008 

-113.5 

0.0003 

-97.8 

0.0001 

-125.5 

1.35565 

1.35565 

110.000 

0.2115 

0.0048 

-100.1 

0.0036 

-133.3 

0.0019 

-121.0 

1.12404 

1.12404 

160.000 

0.3077 

0.0042 

-11.5 

0.0265 

-37.7 

0.0152 

-15.2 

0.98114 

0.98114 

210.000 

0.4038 

0.0164 

-90.5 

0.0536 

-112.6 

0.0610 

-135.7 

0.88177 

0.88177 

260.000 

0.5000 

0.0152 

19.6 

0.1493 

136.0 

0.1379 

-150.8 

0.80754 

0.80754 

310.000 

0.5962 

0.0217 

32.4 

0.2701 

152.4 

0.0923 

174.1 

0.74936 

0.74936 

360.000 

0.6923 

0.0339 

-74.7 

0.1904 

166.9 

0.1795 

112.7 

0.70218 

0.70218 

410.000 

0.7885 

0.1055 

-88.6 

0.0507 

-175.2 

0.3214 

102.0 

0.66292 

0.66292 

460.000 

0.8846 

0.1852 

-90.0 

0.0912 

-15.1 

0.4431 

98.9 

0.62959 

0.62959 

510.000 

0.9808 

0.2619 

-90.2 

0.2137 

-8.4 

0.5413 

97.5 

0.60083 

0.60083 

560.000 

1.0769 

0.3317 

-90.1 

0.3178 

-5.9 

0.6202 

96.7 

0.57568 

0.57568 

610.000 

1.1731 

0.3935 

-90.1 

0.4053 

-4.6 

0.6833 

96.1 

0.55344 

0.55344 

660.000 

1.2692 

0.4477 

-90.0 

0.4787 

-3.7 

0.7343 

95.6 

0.53360 

0.53360 

710.000 

1.3654 

0.4950 

-90.0 

0.5403 

-3.1 

0.7760 

95.2 

0.51575 

0.51575 

760.000 

1.4615 

0.5362 

-90.0 

0.5922 

-2.6 

0.8102 

94.9 

0.49957 

0.49957 

810.000 

1.5577 

0.5721 

-90.0 

0.6364 

-2.3 

0.8384 

94.6 

0.48483 

0.48483 

860.000 

1.6538 

0.6035 

-90.0 

0.6740 

-2.0 

0.8620 

94.4 

0.47133 

0.47133 

910.000 

1.7500 

0.6312 

-90.1 

0.7063 

-1.8 

0.8819 

94.1 

0.45889 

0.45889 

960.000 

1.8462 

0.6556 

-90.1 

0.7343 

-1.6 

0.8987 

93.9 

0.44739 

0.44739 

1010.000 

1.9423 

0.6771 

-90.1 

0.7585 

-1.4 

0.9130 

93.8 

0.43671 

0.43671 

1060.000 

2.0385 

0.6963 

-90.2 

0.7797 

-1.3 

0.9253 

93.6 

0.42676 

0.42676 

1110.000 

2.1346 

0.7134 

-90.3 

0.7983 

-1.2 

0.9358 

93.4 

0.41746 

0.41746 

1160.000 

2.2308 

0.7287 

-90.3 

0.8147 

-1.1 

0.9450 

93.3 

0.40874 

0.40874 

1210.000 

2.3269 

0.7425 

-90.4 

0.8292 

-1.0 

0.9530 

93.1 

0.40055 

0.40055 

1260.000 

2.4231 

0.7549 

-90.5 

0.8421 

-0.9 

0.9599 

93.0 

0.39283 

0.39283 

1310.000 

2.5192 

0.7662 

-90.5 

0.8536 

-0.9 

0.9661 

92.9 

0.38554 

0.38554 

1360.000 

2.6154 

0.7764 

-90.6 

0.8640 

-0.8 

0.9715 

92.8 

0.37865 

0.37865 

1410.000 

2.7115 

0.7857 

-90.6 

0.8733 

-0.8 

0.9763 

92.7 

0.37211 

0.37211 

1460.000 

2.8077 

0.7942 

-90.7 

0.8817 

-0.7 

0.9805 

92.6 


IContainer ship model 
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0 SPEED 
+ 

+ 


0.0000 


WAVE ANGLE 


15.00 DEG. 


VERTICAL PLANE RESPONSES 
(NON-DIMENSIONAL) 


WAVE 

ENCOUNTER WAVE 

WAVE/SHIP 










S U 

R G E 

H E 

A V E 

P I 

T C H 

FREQ 

U E N C I E 

S LENGTH 

LENGTH 











AMPL. 

PHASE 

AMPL. 

PHASE 

AMPL. 

PHASE 

4.49617 

4.49617 

10.000 

0.0192 

0.0000 

-125.6 

0.0000 

-110.8 

0.0000 

-107.4 

1.83555 

1.83555 

60.000 

0.1154 

0.0004 

-176.0 

0.0004 

146.0 

0.0002 

-172.3 

1.35565 

1.35565 

110.000 

0.2115 

0.0037 

-118.6 

0.0045 

-168.1 

0.0020 

-153.9 

1.12404 

1.12404 

160.000 

0.3077 

0.0069 
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